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Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for the IBM Data Processing Systems. If such announce- 
ment indicates a change to the program decks or tapes, a 
complete new program, if needed, should be requested 
from the Program Distribution Center, 



A cknowle dge men t 



This program is a modification of one for the basic 
i620 with paper tape input and output written by 
C* R. Nichols, IBM, Inglewood, calif ornia. The only 
changes made were those necessary to enable his 
program to be used with a card lo20* Much, of the 
information in the writeup is taken directly from 
his writeup* 



TABLE OP CONTENTS 



lm Program description 

A 9 Machine requirements 
Br Error stops 
II* The program 

A« Operating instruct ions 

To load initial data 

To effect cost changes 

To effect requirement changes 

To solve and obtain output 
Bm Data preparation 

Matrix data 

Cost changes 

Requirement changes 
C* Interpretation of results 

Cost or requirement changes 

Dual or simplex algorithms 

Final basis output 

Pinal non^basis output 

Matrix punchout 
Dm Functions of individual subprograms 

Data load routine 

Cost change routine 

Requirement change routine 

Shadow price calculator 



Dual algorithm 11 

Simplex algorithm 12 

Pinal basis output 12 

Non-basis output 12 

Matrix punch routine 12 

E# Use of storage 12 

?• Makeup of program decks 13 

(*• Storage of matrix in memory lij. 

III. sample problem 15 

A* Matrix 15 

B* Matrix in storage 15 

C# Input 1^ 

!• Fixed point 15 

2. Floating point 15 

D» Output 16 

!• Final matrix punched 16 

2. Typewriter output 16 

IV* Method of loading program 16 

Mm Flow charts 17 

km General flow chart 17 

Bm Data loader 19 

C# Cost changer 22 

D* Requirement changer 2l|. 

E# shadow price calculator 27 

F^ Dual algorithm 31 

G. Simplex algorithm 35 

H* Final basis output 39 



I# Fix subroutine 

JY Non-basis output 

Km Matrix punch routine 
VI • Listing of program decks 

A* Data loader 

B* Cost changer 

C# Requirement change r 

D# Solution deck 



I 



T# PROGRAM description 



This program solves linear programming problems with 
output of detailed re suit s* That is, given coeffi- 
cients &±i$ cost coefficients c-j, and requirements 
b± 9 determine such that 

^ a ij x j * b i with xj^O 
j 

and 

£1 c jxj » maximum 
J 

Computations are performed by the Dual Algorithm until 
a feasible solution is obtained* Control is then given 
to the simplex Algorithm for optimization* Cost changes 
and requirement changes can be made after loading 
original matrix or after solving original matrix* 



JU Machine requirements 

1* Basic 1620 with card Input and output* 

2* Any size storage can be used* The larger the 
storage, the larger the problem that can be 
solved* The size of the problem is restricted 
by the following relationship: 

(m + 2) (n + 3) ^ memory - 3920 

*TTT 

where- m Is the number of restrictions 

n is the number of no n-ba sis independent 

variables 
memory is 20,000, lj.0,000, or 60,000 

3* Source language: actual machine language 

i}_* subroutines: floating point subroutines supplied 
with program* No other subroutines needed* 

Input media: card reader and console typewriter* 

6* Output media: card punch and console typewriter* 

7* Running time: The precise time required per 
iteration depends on the size and density of 
the matrix* As an approximation, a problem 
with 30 equations and I4.0 non-basis variables 
requires about 20 seconds per iteration* 

8* Decimal accuracy: All computations are performed 
In 2-and-8 floating point form* Matrix input 
can be either fixed point or floating point* 



Error stops 



1* All stops may be recognized by displaying 
IR-1 in MAR* MAR will contain the address 
of the halt instruction plus 12* 

Halt instr* Reas o n 

01126 Floating point overflow 

01288 Floating point attempted division 

by zero 

00012 Loader - ready for next code 

02I4.58 Cost changer - Change entered for 

a non-existent identification 

number 

00012 Cost changer - ready for next code 

O31IO RHS changer - change entered for 

a non-existent identification 

number 

00012 RHS changer - ready for next code 

02E>3t|. Dual algorithm - inconsistent matrix 

026£k simplex algorithm - unbounded solution 

02838 Matrix punch - ready for next code 

2* Error procedure 

a* cost changer ( 02ij.£8 ) • Erroneous identification 
number is the last one typed* 

If data are being entered through the type- 
writer, depress RESET and INSERT* 
Type ij.9 02126 • 
Depress RELEASE and START* 
Corrected or new change data may now be 
entered normally* 

If data are being entered from cards, the 
erroneous datum may be ignored by following 
the same steps as above (for entry through 
typewriter )• 

If it is desired to enter the change correctly 

through the typewriter, the following steps 

must be taken: 

Set console switch two on* 

Deoress RESET and INSERT* 

Type lj.9 02126* 

Depress RELEASE and START* 

Enter correct datum through typewriter 

(complete with record mark)* 

Depress RELEASE and START* 

Mien typewriter is again selected, depress 

RESET. 

Set console switch two off* 



Depress INSERT, VZ^?ASE 9 and START. 
Processing will n v continue normally, 
reading data from cards • 

b. RH3 Changer (03110): Erroneous identification 
number Is the last one typed * 

If data are being entered through the type- 

x;riter, depress RESET and INSERT* 

Type 1l9 02126 s 

Depress RELEASE and START. 

Corrected or new change data may now be 

entered normally. 

If data are being entered from cards, the 
erroneous data may bo ignored by following 
the same steps as above (for entry through 
typewriter). 

If It is desired to enter the change correctly 
through the typewriter, the following steps 
must be taken: 

Set console switch three on. 
Depress RESET and INSERT ♦ 
Tyne l< r 9 02126. 
Depress RELEASE and START. 

Enter correct Identification number digits 

plus record mark) through typewriter. 

Depress REIJSASE and START • 

Enter requirement change (including record 

mark) through typewriter. 

Depress RELEASE and START* 

When typewriter Is again selected, depress 

reset/ 

Set console switch three off. 
Depress INSERT , RELEASE, and START* 
Processing will now continue normally, reading 
data from "cards. 

c. Dual algorithm (0253U-) s To obtain solution 
existing at the time of discovery of inconsis- 
tency: 

Run otit cards from reader. 

Replace cards In reader hopper starting with 

Pinal Basis Output program. 

Depress READER START. 

Depress RESET and INSERT. 

Type l|_9 0002k. 

Depress RELEASE and START* 

d. simplex Algorithm (026£if.): To obtain solution 
existing at the time of discovery of unbounded 
solution: 



Depress RESET and INSERT* 

Type l-[.9 0002k. 

Depress RELEASE and START* 

e# Overflow-underflow: Overflow causes a 

stop at 01126 or 01288 . Underflow causes 
the result to be set to zero, and process- 
ing continues* Either condition indicates 
6i$her machine malfunction or unreasonable 
data* 

f • Any other error: unrecoverable* 
Restart problem* 



THE PROGRAM 



Operation Instructions 

1* set typewriter margins at 12 and 92; tabs at 
2ij., 36, k Q 9 60, and 69* 

2* It is important to note that no more than one 
instruction at a time may be inserted at any 
time after the initiation of the Loader routine* 

3* To load initial data: 

a* Ready the LOADER deck in the card reader 
(including matrix input data)* Depress 
READER START* 

b* Depress RESET and INSERT* 

c* Type 16 00010 00000. 

d. Depress RELEASE and START. 

r $iis initiates memory clear' routine* 

e* After at least one second* depress INSTANT 
STOP, RESET, and INSERT* 

f* Type 36 03826 00500 
1;9 03826 

g* Depress RELEASE and START* 

h* Program deck loads, data are loaded, and a 
halt at location OO03.2 is executed. Opera- 
tion may now proceed to the COST CHANGER, 
RHS CHANGER, or SOLUTION programs • 
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To effect cost changes: 

a. Ready the COST CHANGER deck in the card 
reader (including data change cards if 
changes will be made from cards)* Depress 
READER START . 

b* set console switch two off if cost changes 
will be entered from carcTreader; switch 
two on if data will be entered through the 
typewriter* 

c* Depress START* 

d* If data are entered through the typewriter, 
RELEASE and START must be depressed after 
each entry. (Don*t forget the record mark 
to terminate each entry*) 

e* After all data are entered a halt at location 
00012 is executed* Operation may now proceed 
to the RHS CHANGER or SOLUTION programs. 

£. To effect requirement changes: 

a* Ready the RHS CHANGER deck in the card 
reader (including data change cards if 
changes will be made from cards)* Depress 
READER START. 

b* set console switch three off if requirement 
changes will be entered from card reader; 
switch three on if data will be entered 
through the typewriter* 

c* Depress START. 

d. If data are entered through the typewriter, 
RELEASE and START must be depressed after 
the entry of each identification and each 
change. (Don't forget the record mark 
to terminate each entry. ) 

e* After all data are entered a halt at location 
00012 is executed* Operation may now proceed 
to the COST CHANGER or SOLUTION programs. 

6* To solve and obtain output: 

a. Ready SOLVER deck in the card reader* 
Depress READER START* 



b. set console switch one on to type iterations, 
off to omit iteration typeout* 
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c* set console switch four on to obtain punch- 
out of final matrix, offT o suppress punch- 
ing* 

d. Depress START. 

e. After all output is complete, a halt at 
location O2838 will be executed* Operation 
may now proceed to the LOADER, COST CHANGER, 
or RHS CHANGER programs. 



B* Data Preparation 

1* The first card to be entered is a c ase identification 
card* A five digit alphabetic identification ~~ 
is punched in columns 1 through £* columns 6 
through 80 are blank* 



2v The second card specifies the size of the problem 
being entered* It is punched as follows: 

Information Card columns 



(zero with a minus sign) 1 
number of equations (3 digits) 2 - I}. 

(zero) 5 
(zero with a minus sign) 6 
number of non-basis variables (3 dig* ) 7 - 9 
(zero) 10 
input code 11 
a* 1 for row- column, fixed point 
input 

b # for a complete floating point 
matrix 

Columns 12 through 80 are blank* 

The matrix can be entered in one of two forms? 
floating point (the entire matrix is entered* 
column by column, exactly as it is to be stored* 
including zero entries ) or fixed point (entries 
can be in any order, accompanied by their res- 
pective row-column designations* and zero ele- 
ments need not be entered)* 

For floating point entry, eight elements are 

punched on each card, with a minus sign over 

the first position of each ten position field* 

and a minus sign over the last position of 

negative fields (units position)* There are ten 

positions per field, with the first two positions 

specifying placement of the decimal point (5>0 . |^ ) 

designates a decimal point to the left of the 
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first non-zero digit; £l designates one non- 
zero digit to the left of the decimal; i{.9 desig- 
nates that the first non-zero digit is in the 
hundredths position, etc.) and the last eight 
digits (with no leading zeroes) giving the 
element itself* A zero element is specified 
by ten zeroes, with a minus sign over the firsts 

The elements must be punched in column order, 
starting with the first basis cost. See the 
matrix layout included with this writeup to 
know how the elements are stored in the matrix. 

£♦ Fixed point matrix elements are entered four 
to a card. 

Two types of cards are necessary for matrix 
elements : 

a. Rox-j- column cards designate where the elements 
will be stored in the matrix. Row-column 
designations are punched- in columns 1-6, 
21-26, ltl-^6, and 61-66 in the form RRRCCC 

(a minus sign accompanies the first digit 
of each row and column designation). The 
column designation specifies column number 
considering only non-basis vectors. The 
requirement variables have only a row desig- 
nation, and 000 should be entered for their 
column. 

b. Matrix element cards contain fixed-point 
entries including optional leading sign, 
up to eight numeric digits, and mandatory 
decimal point. Each entry must be terminated 
by a record mark (0-2-8 punch). The four 
elements are punched beginning in columns 

1, 21, );.l, and 6l for each card. 

Elements may be entered in any order, but corres- 
ponding row-column designations and matrix 
entries must be in the same positions of suc- 
cessive cards (row-column card preceding). 

Loading of zero elements is optional 



If any card contains fewer than four elements, 
the first blank field (column 21, i|.l, or 6l) 
of the row-column card must be -punched with a 
record mark (0-2-8 punch). 
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Item Description Row-column 



aj[j matrix element of the Ith row . 
and jth non-basis vector xxxxxx 

*■ 

h± element in the ith row of 

the requirement vector xxx$00 

6* Fixed point cost entries are also four per card, 
but the row-column designation is on the same 
card as the cost* 

Entries consist of a six digit row- column desig- 
nation followed by a ten digit ID/cost, where 
the ID is a four digit numerical designation 
to identify the cost variable, and cost is a 
six digit fixed point entry with the decimal 
after the third digit* 

Item Description Row~ column 

ci cost per unit of basis — % 

variable for row i xxx^OTxxxxxxxxxx 

cj cost per unit of jth _ .^a-^-*-^, — ^ 

non-basis vector OOlxxxxxxxxxxxxx 

negative cost is signified by a minus sign over 
the units position of the cost field* 

If any card contains fewer than four elements, 
the first blank field (column 21, I4.I, or 61) 
of that card must be punched with a record mark 
(0-2-8 punch)* 

7* The^last data entry for fixed point mode torus t 
be 000 to terminate the loading operation* oSiis 
can be punched in columns 1-3 of a new card, 
or in 21-23, ijl-ijj, or 6I-63 or the last row- 
column or cost entry card* 

8* Entry of cost changes: 

Data for the COST CHANGER routine can be entered 
either from cards or from the typewriter* If 
entered from cards, one entry is punched in 
each card, in columns 1-10* The entry consists 
of the four digit ID and six digit cost figure , 
with a minus sign over the first digit of the * 
cost (column 5) and over the last digit of negative 
costs (column 10)* Following the last cost 

change card must be a card punched with four zeroes £>) 
(columns 1-lj.) to terminate the operation* 1> 
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If entered from the typewriter, the format. Is 
the same as from cards: four digit ID and six 
digit cost figure with flag over the first digit 
of all cost figures, and over the units position 
of negative costs • The last entry must be four 
Eeroes to terminate the operation* 

9* Entry of requirement changes: 

•phere are two entries for each requirement change: 
the identification and the chan g e in the value 
of that requirement element. Wte that the 
now requirement element is not entered, but 
rather the difference between the old and new 
requirement elements* If changes are entered 
from cards, the identification and change are 
entered on separate cards, one per card* 

The identification is entered first, and con- 
sists of the ID associated with the cost element 
for that row* Thus, to change bi ■ b3, use the 
ID associated with c^ - co* No flags or record 
marks are needed in entering the ID either 
from cards or typewriter* It is punched in 
columns of card entries, with £-80 blank* 

After entering the ID, its associated requirement 
change is entered (by typewriter or on a separate 
card)* The entry starts in column 1, if a 
card is being used, and consists of optional 
leading sign, up to eight numeric digits, and 
mandatory decimal point* It must be terminated 
by a record mark, whether it is entered from 
cards or from the typewriter* 

The last entry must be four zeroes to terminate 
the operation* 

C* Interpretation of Results 

1# Cost changes or requirement changes are typed 

if the input is from cards* If input is through 
the typewriter, the input operation itself 
produces a log of the changes made* 

2* Dual Algorithm and simplex Algorithm: 

Results of iterations are typed only if console 
switch one is on* If it is desired to monitor 
the course of the solution onlh occasionally, 
this may be done by setting switch one on and 
off periodically* The following points should 
be borne in mind; 
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a* Page headings for iteration output are 

typed only if switch one is on at the be- 
ginning of the solution* 

b* The iteration count is reset x*hen the dual 
algorithm is finished and control passes 
to the simplex algorithm* 

[Che information which may be typed includes 
the iteration count, the current value of the 
profit function (floating point), the last 
variable to leave the basis, and the variable 
which entered the basis* 

Pinal Basis Output: 

The final value of the profit function is given 
in floating point* 

For each final basis variable, the following 
information is supplied (in fixed point): 

a. Identification/cost coefficient* 

b* Activity* 

c* The limits of the cost coefficient over 
which the current solution is optimal* 

d* The variables which limit the range of the 
cost coefficient and will enter the basis 
if a limit is exceeded* 

Final Non- Basis Output: 

For each non-basis variable, the follox^ing 
information is given (in fixed point): 

a* Identification/cost coefficient* 

b* shadow price - the penalty to the total 

system if a unit of this variable is forced 
into the final solution* 

c* The limits of activity over tirhich the shadow 
price applies* 

d* The variables which limit the range of 
applicability of the shadow price* The 
upper limiting variable is the one which 
would leave the basis if the associated 
non-bftsis variable were forced into the 
solution* 
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% A punchout of the complete final matrix is op- 
tional under control of console switch four* 
The cards punched also contain all necessary 
supplementary information and are in a format 
which is suitable for direct re-loading by 
the LOADER routine. 

Functions of Individual sub-programs 

1. The data load routine reads and stores parameters 
for the problem. It then proceeds to load and 
store the matrix in the proper format. In the 
case of fixed point input* storage locations 

are computed from the row/column designations $ 
and the input elements are converted to floating 
point notation as they are read and stored. 

2. The cost change routine reads ID/cost elements* 
searches the matrix for a corresponding iden- 
tification number, and stores the new cost 
coefficient in the proper location. 

3* The requirement change routine reads variable 
identification numbers and changes for the 
associated requirement element. The change is 
converted to floating point notations, the 
matrix is searched for the identification num- 
ber, and the requirement vector is updated. 

k. The shadow price calculator converts all cost 
coefficients to floating point notation and 
evaluates Zj-cj for each non-basis variable. 

5>. TJie dual algorithm computes a feasible solution 
for the problem. Its operation may be monitored 
allowing the typing of the value of the functional 
and the identifications of the variables that 
enter into and are removed from the basis. It 
should be noted in this regard that the objectives 
of feasibility and optimality may not be compati- 
ble at this point, and therefore the functional 
may actually decrease from one iteration to the 
next while the dual algorithm is in control. 

During the reduction of the matrix of coefficients 
a test is made to determine whether the resulting 
quantity has a magnitude less than an amount 
called "essential zero." Any element falling 
into this category is set to zero in order to 
avoid computation with numbers which are actually 
not significant. "Essential zero" has been set 
to 10"°. Its value (the corresponding floating 
point exponent) is located in memory addresses 
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31^ and 31^5 and may be changed if desired* If 
this change is made, it is important to note 
that the high order digit must carry a flag*. 

6 # The simplex algorithm starts with any feasible 
solution and computes the optimal feasible solu- 
tion. Its operation may be monitored by allow- 
ing iteration typeout as in the case of the dual 
algorithm. 

iJxe discussion of essential zero as applied to 
the dual algorithm in f>. above is equally valid 
for the simplex algorithm. In this case the test 
value is located in memory addresses 326I4. and 
326£. 

In the case of a tie in the quotient used as the 
criterion for choosing the variable to leave the 
basis on any iteration, the choice is resolved 
by using the largest denominator as a secondary 
criterion. 

?• The basis output routine searches out and com- 
putes the information described in section C# 
above. Floating point quantities (with the ex- 
ception of the functional) are converted to 
fixed point before being typed. 

8. Hie non-basis output routine searches out and 
computes the information described in section 
Cm above. Floating point quantities are converted 
to fixed point before being typed. 

9# The matrix punch routine resets the permanent 

instructions in low memory to allow re -initiation 
of the input or change phases. Problem parameters 
and the final matrix are punched into paper tape 
if called for by the setting of console switch 
four. 

E. Use of Storage 

Memory locations from 00012 through 03919 are used 
by the program. Memory from 03920 up to the re- 
quired location is used for the matrix and its mani- 
pulation. Locations 00000 through 00011 are avail- 
able for the insertion of one instruction at a time 
through the typewriter. A halt instruction in posi^ 
tion 00012 ensures that only this one instruction 
will be executed. 

00012-000lt2 load routine 
OOOL0-OOO58 case identification 
00060-00061}. S address 
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0006^-00069 M and 2 

00070-000711. N 

00075-00079 M 

OOO8O-OOO99 product area 

00100-00l}.00 arithmetic tables 

00li.02-0206Ii floating point subroutines 

variable programs 

O3758-O3917 input area (record mark in 0391o ) 

03920- matrix 

Program Decks 

1* LOADER 

3 r>re-load cards* 
Arithmetic tables* 
floating point routines* 
Data load routine* 
Matrix data to be entered. 

2. COST CHANGER 



2 pre-load cards* 
Arithmetic tables. 
Floating point routines* 
Cost changer routine* 

Cost change data (if entered from cards). 
3* RHS CTIA?!GER 

2 pre-load cards. 
Arithmetic tables. 
Floating point routines* 
Requirement change routine* 

Requirement change data (if entered from cards > 
ij.. SOLUTION 

2 pre-load cards* 

Shadow price calculation routine* 

Page headings* 

2 pre-load cards* 

Dual algorithm routine* 

2 pre-load cards* 

Simplex algorithm* 

2 pre-load cards* 

Final basis output routine. 

2 pre-load cards* 

Fix and print routine* 

Page headings* 

2 pre-load cards* 

Non-basis output routine. 

Page headings. 

2 Pre-load cards. 

Matrix punch routine* 



G# storage of matrix in memory 

The matrix is stored in column sequence and uses 
memory from 3920 upward to the extent required by 
the problem size* The following diagram indicates 
the manner in which matrix elements are stored: 



3920 






S+10(2m+3) ID/Cj S+10[(n+l)(mf2)-lj 




S 


Functional 


S410(2m+4) S4loj<n4l)(m42j] 


Working 
column 

■ 

3920+ 
10(m+l) 


ID/d 
S-llOm 


S410(m-f3) 
S410(2m+2) 


S410(2m+5 ) a+iofdifl ) (m+2 

Increasing j 

a ij: S4l0[(3tl)(mf2)ti] 

Incpeasing 
i 

■ 

S+10(3m44) Sfl0£(nf2) {n f 2)-2] 



S is computed for each problem and is equal to 3920 
4 10(mf3)« All addresses are field addresses* 
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o 



III. SAMPLE PROBLEM 



A. Matrix 



Requirement 

Basis variables Non-basis variables vector 



9901 


9992 


0001 


0002 


9902 




1 




1 


2 




6 




I 


1 


-1 


-1 i 


-k 





-10 


1 


1 


° ! 





Variable 
ID 



Cost 



B* Matrix in storage 



O 





0000000000 


"5001001000 0*002001000 9902000000 




0*000000000 


TTOOOOOOOOO "0*000000000 0*000000000 


•5901000000 
9992010000" 


51 60000000 

511^.0000000" 


"fllOOOOOOO ">120000000 0000000000 
5110000000 "TllOOOOOOO "fllOOOOOOO 



2X3 
S00200"00301 
001001 

l.fc 

002(503 

750100T2901000000 
*50IO"039902000000 

0*00 



input 

!• Fixed point 



0*010*02 
2.* 
0*010*00 
6.4 

T|02"50"*5 , 992 01000*5 



0*020*01 
i">02"500 

"5oTooiT5ooiooiooo 



0*020*02 

1.* 

0"oTo"02*5o02001000 



2. Floating point 



2X3 



0000*000300 

" aQDQH 99g0l0nnnnn D nnnnnnnnn0n00 000N160000o00N 

nL^ nn T !!!!nL « onnooo2nolooonooonooooON,2oooonnnN ^ 



r 
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D# Output 

1# Pinal matrix punched 



2X3 

0002.0000300 

R902000000 0001 00 1 OOOOOOOOOOOOON 1 60O00000M21 ooooooon I6ooooononorv2 OOlOOONl ioooaooo 
N 1 30000000N 1 20000000R9920 1 0000N2 1 0O000O0N 1 1 0000O0000000000O0R9O 1 OOOOOON 1 1 0000000 
Nl 1 ooonooONl 1 oooooooo oooooo oo nooooooooooooooonnonnnnooooonnnnnnooooooonnnnnnon 

2# Typewriter output 



CASE 



2X3 



ITER NO 


FUNCTIONAL 


501 


51 0000000 


502 


5130000000 


501 


5160000000 


FUNCTIONAL 


5160000000 


VAR/COST 


ACTIVITY 


9902000000 


10.0000 


0001001000 


6.0000 


VAR/COST 


SHAD PRICE 


5002001000 


1.0000 


§992010005 


10.0000 


9901 000000 


1.0000 



VAR OUT 

799201 0005 
^901 000000 
1002001000 



LIM VAR 

5002 
5002 

LIM VAR 

5000 
§902 
0000 



VAR IN 

5002001000 
§902000000 
0001001000 



LOWER LIM LIM VAR UPPER LIM 



.3333- 9992 
.5000 0000 



9999.9000- 5001 
10.0000- 5000 
9999.9000- 5001 



10.0000 
9999.9000 



LOWER LIM LIM VAR UPPER LIM 



3.0000 
9999.9000 
6.0000 



IV. METHOD OP LOADING PROGRAM 

The "LOAD" key cannot be used to load programs in this 
routine since the first 80 positions' of memory are not 
free. Therefore cards are read into locations 03838 
through 03917* with a record mark in 03918 being utilized 
to transfer up to six instructions to their place in 
memory. The first two instructions on each card con- 
tain the "transfer record 11 instruction and a branch to 
location 03826, which contains the read instruction 
36 03838 005"00. The first two cards of each routine 
perform the function of setting th» read instruction 
in 03826 and the record mark in 03918. The first card 
is read into 03826, and control is transferred there, 
using Instructions 36 03826 00500 ij.9 03826. This 
must be done by the operator a* the beginning of the 
"LOADER" program, it is done automatically for all 
following routines* 



COMPUTER 
TECHNOLOGY 
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GENERAL PLOW 
page 1 



MANUALLY 
CLEAR 
MEMORY 

I 



READ IN 

ARITH. 
TABLES 



PAD 

FLOATING PT 
SUBROUTINES 



iHEAI 
LOAI 
PROC 



ID IN 
I LOADER 
PROGRAM 



I 



(re 

MA 1 



[READ IN 
TRIX DATA 



FLOAT & 
STORE IN 
MATRIX 



COST CHA NGES, 



RHS CHANGER 



READ IN 

ARITH. 

TABTiKg 



READ IN 

FLOATING P 
STTBROTTTTTJR 



1. 



READ IN 
COST CHANGE 
PROG_BAE 




READ NEW 
COST FRO* 
TYPEWRITER 



READ NEW 
COST FROM 

£AJQa_ 




RHS 



CHANGER 



SHADOW PRICE 



SODUTION 



(READ IN 
SHADOW PRICE 

PR OGRAM 
T 




f READ 


IN 


PAGE 




TTRAU* 


IT?GS ... 





J TYPE 
PAGE 
j HEADINGS 



I 




fREAD IN 
DUAL ALSO. 
PROGRAM 



TE 

FEASIBLE I 
SOLUTION 




READ IN 

ARITH. 

TABLES 



floJJtBIg pti. 
subroutines 



I 



LD IN 
*HS CHANGE 
PROGRAM 




READ NEW 
REQ. FROfe 
TYPEWRIT !R 



L 



PLACE NEV 
REQ. IN 
MATRIX 

7~ 



READ 2S W „ 
REft. FROM 

CARDS 

_ 



\ 



SHADOW PRICE 



CoIt CHANGER 
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GENERAL PLOW 
Page 2 



1 



[READ IN 
SIMPLEX ALCO. 
JEBQGBA1I J 



I COMPUTE 

OPTIMUM 
ITSASIBLE 
HflLTTTTfiN 



s 



HEAD IN 

BASIS O/P 
PROGRAM 



READ" FIX 
SUBROUTINE 



I 



f READ IN 
PAGE 



TYPE 
PAGE 
HEADINGS! 



TYPE 
BASIS 
OUTPUT 



I 



READ IN 
KON-BASIS 0, 
PROGRAM 



f READ IN 
PAGE 

HEADINGS 



FT 



PAGE 
HEADINGS! 



NON-BASIB 
OUTPUT 



READ IN 
MATRIX PCH. 
PROGRAM 



PUNCH 
«ATRIX 
SW 1; ON 



3 



COST CHANGER 



LOADER 



RHS CHANGER 



r 



2259 

flNITIALlfeE 
STORE 
INSTR. 



V 2270 

{READ PL. 
! INPUT 
[ CARD 

f 2282 

STORE 

IN 
MATRIX 

t 22% 
MODIFY 
PREV. 
INSTR. 



~ J 

,NO Y-^?* 
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CARD 



— (YEEJ 




DATA LOADER 
Page 1 



_S3L62 



READ M, 
N, AND 
INPUT 



CODE 



Z1..23L' 
ICALC. & 
I STORE S» 
*f*10(M+3)j 

STORE I 
M, M+2, 
N 




t 22U.6 



CODE=0~ 



ho) 



READ 
ROW/COL 
CARD 



2326 




/ 2338 



^ A2338 
P ^YE9-^ROW=^ >-(NO 



COMPUTE 
MATRIX SIZE 

10(N42)(M42) 



5362. 



COMPUTE & STORS 
STORAGE LOCATION 
S^10t(j+1)(M«-2H] 



SET FIELD 
FLAGS FOR 
ELEMENTS 




W-^>-^Es)- 



V 12 



HALT 



I 



READ A 
MATRIX 
CARD 



TRP 1st I 
COST 



READ 
NEXT 
CODE 




g. 3 
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DATA. LOADER 
Page 2 



2$22 



SET 
i = 1 



LafifaA 
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DATA LOADER 
Page 3 



292$ 
SET £«1 

SET K=9 



*3- 



, 

XHITIALIZI 




INSTR. 



J 



297^ _ 
2»K 

, 

'MOVE Dl^kl 1 
TO PLOIWAREA 



K = K-l 




$isw k=o? 



3086 



i 313Q 



SET 
PROPER 
SIGN 

3l9i 



V 3098 

SET 
DIGIT i 
TO ZERO 




■O 1 



STORE 
ELEMENT 
IN MATRIX 



311,0 

1 = 14-1 



T 





I SLIDE 
* READ-IN 
AREA 




RESET 25^0 
TO READ 
IN STR. 

r" 



O-OP REAj 
INSTR. 
L2SLQ1 




Pg.l 



( 2338 \ n 
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COST CHAN3-ER 
Page 1 



HALT 



o 



2126 



I PR 



. ^\21 
JWITC§- 
,'NO 2 ON ? 

j 2 138 




2126 



-/EES 



D ID/ COST 
FROM CARD 



2k 
READ 
NEXT | 
CODE J 



JjL217Q 



READ ID/COST 

FROM 
1 TYPEWRJTgR.,.. 




TYPE 
ID/COST 



V2266 



SET K2M 
SET i«l 
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COST C SANGER 
Page 2 




SET K=N 
SET j*l 




2USP 



SET 
ADDRESS 
OP ID^ 




T2H92 



STORE 
NEW 
COST 



6. 1 



SET 
ADDRESS 
OP IDj 
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RHS CHAJDER 
Page 1 



READ' ID 




FROM 




CARD 


1 1 




COST 
CHANGER 




/On. 



o 
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© 



2h< 



SET 
i - 1 



RHS CHANGER 
Page 2 




EXP«EXP+t 



Jt22Q2 



SET k»9 
SET i»l 




INITIALISE 
INSTR. 



MOVE DIG| 
i TO 
fLQ AT AREK 



k s k-3 



2314. 



i = 1+1 
k = k-1 








SET 
DIGIT 1 
TO ZERO 



2906 



P set 
; PROPER 
: SIGN 



8 i+1 



T 



2930 

SET k*ui 
SET i"l 



T" 
© 



SHIFT OF 
DECIMAL 
POINT 




□FT OOT 
(LEADING I 
! ZERO I 

t267Q 




RES CHANGER 
Page 3 



HALT 



^NO )— c<0 ?^-{YE2| 




o 



i 



221,0 
SET 
WORK AREA 
TO ZERO 



^0 
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SET i«0 
SET x-m 



L2150 



SHADOW PRICE CALC, 
Page 1 



i ■ 1*1 



j 2162 



MOVE Ci 
TO 

UORK AREA 




SAVE 

SIGN 




dXP"EXP-l 



SAVE k 



f 

Formalize! 



PILL OUT 
TRAILING 
ZEROES 



SET 



1 







F~STOHE I 
FLOATED 0{ 
g| WORK C^ L. 




SET 
WORK AREA] 
TO ZERO 



© 
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SHADOW PRICE CALC. 
Page 2 



x - x-l 




MODIFY LOOP 
TO FLOAT & 
NEGATE Cj'3 



SET J«0. 
SET x«n 



E 



2150 




32162 



MOVE 04 
TO J 

Work areai 




2230 



SET 

EXP = 53 
k = 1 



t229Q 



SAVE 
SIGN 




FILL OUT 
TRAILING 
ZEROES 




g. 1 



o 



o 
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SHADOW PRICE CALC. 
Page 3 



V 2lt66 



STORE FLOATED 
3 4 IS SHADOW 
PRICE ROW 

_,. .^2fr.90 



x = x-1 



~A?502 
(®) — <x ■ 0? 




g. 2 



SET x=n 



SET 
FUNCTIONAL 
TO ZERO 

/ T - 

ISET 3*0 



3 



SET 1=1 
SET y=m 




2752 





-30 



SHADOW PRICE CALC. 
Page k 







J2212 



X ■ x-1 




READ 
PAGE 
HEADING* 




DUAL 
ALGORITHM 



;pg. 
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© 



© 



© 




2iik 



SET baia 


» 


SET z ■ 


m 


SET 1 ■ 






5*2 



1 « 141 




x ■ x-1 




SE £ R max 

S -09 

SET y*»n 
SET 3=0 



SIMPLEX 
ALGORITHM 



SET ROW t 
ADDRESS 



12330 



3 « J«. 



6 



DUAL ALGORITHM 
Page 1 




2330 




(INCONSISTENT 

MATRIX) 



DUAL ALGORITHM 
Page 2 



21150 



3AVE a r j ADDR. 



y - y-i 




SET x*m 
SET 1»0 





DUAL ALGORITHM 
Page 3 



SET j"0 
SET y"n 




SET 1=0 
SET J«0 
SET x*m 




DUAL ALGORITHM 
Page k- 



3GQ2 



SET j«0 
SET ys« 





33SO 



STEP 
tTERATIO! 

coxniT 



1 





g. 1 



KPEj ITERATION, 

<WCTI0NAL, 

ID/cOSTv 

n^osTp 
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S IMP LEX ALGORITHM 
Page 1 



( 2126 J 



91 PA 

SET y ■ n 
SET J > 1 



d mi 



v NOW.^<d m i^fES) 



2210 



SET a^sdj ; 

••"7 J 

H 




SET R mln = <* 
SET x • m 

.SETjLJLO 



2lt 
READ 
NEXT 
CODE 



V 

i = i+1 



25^2 



BASIS 
OUTPUT 



(!) 




CALC 
l/a pk 



SET 
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o 



NO 



SIMPLEX ALGORITHM 
Page 3 



[2&0 



MOVE 



SET ajic 
■ 



I 



X * x-1 



882 




<0? 



2fla 



SET 
a rk « 1 



I 



SET j«0 

set y=n 



a rj" a rj/ a rk 



SET 1-0 
SET j»0 
SET 3C-BI 



/ 





SET j-0 
SET y«n 



^ — - 



a ij"aij-(arj)(aik) 




SWITCH 
ID r & IDj 
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READ 
FIX 
SUBROUTINE 



PINAL BASIS OUTPUT 
Page 1 



TYPE 

yuNCTio: 



ID & 
I TYPE pa. 
IDINCrS 



V 277k 


SET 




SET 





2906 



3ET UPPER 
3ET LOWER 
SET j = 1 
3ET y ■ n 




2906 




SAVE IDj 



TL9 
SE» 1 
iiOWER ■ R 





2202 


SAVE 


3>jl 



«3- 



V321k 



3 = 3*1 
y ■ y-i 



T 




JL33SS. 



UfYPB ACTIVITY 



I 



TYPE LOWER 



ER 33)} 



SET C=C1 



I 



OALC. & TYPE 
LOWER LIMIT 
■ C- (LOWER) 



FIFE UPPER IE 



I 



CALC. & TYPI 
UPPER LIMIT 
■ C- (UPPER) 



I = i-1 



PINAL BASIS OUTPUT 
Pftg« 2 





NON-BASIS 
OUTPUT 



-in- 







FIX SUBROUTINE 
Page 1 



SET EDIT 
WORD IN 
CUTPUT 



ARIA 




SET 
KINUS SlflN 





FIX SUBROUTINE 
P«g« 2 



TURN ON 
BUP FIAO 




pT BLA] 
IN OUTPUT 
P3BIT J 



1 » 


1*1 




\ 




i • J*l 



K2276 
WOVE d jgij 

1 TO OUTPUT 
DIGIT 1 • 





TURN OF 
ALT FLAG 



3 



TAB 
fXPEWSTEE 



TURN OFF 
SUP FLAG 



I 



TYPE 



SET x"5 



OUTPUT 

? 



1 BRA 
1 BA CK 




•*3- 

TYPE 
PAGE 
HEADIHBS 

X 



NON-BASIS OUTPUT 
Page 1 



SET J«l 
SET y«n 



SET UPPER" °* 
SET LOWER*-' 
SET in 
SET x"m 



_^2802^ 




HPS ID/Cj 



I 



TYPE 
SHADOW 
PRICE 

J 



NGN-BASIS OUTPUT 
Page 2 




TYPE 
LOWER ID 

TYPE 
LOWER 
LIMIT 

♦ 

TYPE 
UPPER 1) 

TYPE 
UPPER 
LIMIT 

3 » J+l 

y • y-i 



i 



READ 
NEXT 
CODE 



2t 



MATRIX 
PUNCH 



-W 2802J 



MATRIX PUNCH 
Page 1 



LOADER 




SWITCH 
k ON ? 




I RHS 
^ CHANGER 

COST 
CHANGER 



1 267Q 



PUNCH 
CODE i 




COMPUTE 
MATRIX 
LIMITS 



4.27.78 



i PUNCH 
|L CARD OF 
i MATRIX 




c 



c 



PRE-LOAD* 
CARDS 



ARITHMETIC 
TABLES 



FLOATING 

POINT 

ROUTINES 



360383800500 
1 5039 1 80000Z4903826 
3 1 1 ZQ3 8 a74 9.0 3E2 £l 
31 00 1 00038584903826 
3 100 1 60038584903826 
31 00220038584903826 
.3_la 028 >00 38584 9Q3_826„ 
31 00340038584903826 
31 00402038584903826 



31 00462038584903826 
31 00522038584903826 



4 8 000000000036 Q 30£!600g0 4 Se 3Qga 

0000000000001 0203040002040608000306090210040802161 0050015102 
6CL21_8 1 42QQ7P41 1 PR? 0806 142 2 3 009081 72630000000000 5 060708CT90 
012141618151811 242720242822363520353045403632484455324946536 
0,4846546^7^4453627 1 dp 1 ^4^6789 1 234567890234567690J34 567990 JK 
4567890 JKL567890JKLM67890JKLMN7890-*KLMNO890JKLMNOP90JKLMNOPQ 
L6Q0536QQ53fl49n04340l6005360055flP60l 197qoooo?6q1 1 12QQ44516Q1 



3100582038584903826 
3100642038584903826 



3100702038584903826 
3 1 00762038584903826 
31 00822038584903826 
31 00882038584903826 



3100942038584903826 
3101 002038584903826 



3101 062038584903826 
3101 122038584903826 



3101 182039004903826 
3101 2000385849038 26 
3 1 1 26O038584903826 
31 01 320038584903826 



124000001 101 124000014600494014004901 128004692601 176000002601 
1870182749000000 2201 1 860 1 1 764900570026O1 186Q1 1 761 500847nonn 1 
160085701 187260109201 1892201 1 780 1 1 8947007020 1 1002601 17601 197 
260 1 1 970 1 1 86260 1 1 860 1 1 762 1 1 092 1 1 8926Q 1 1 78n 1 1 S932Q 1 1 78nnnnn 
1501 198000001401 1 780000Q470 1 0940 1 30044007740 1 1 97320 1 1 98OOO0O * 

330 1 1970 0000210085701 I7844008220ii8633nl i86nnonni snnfl47nnnnp 7* 
1501 178000001501 1890000Q2101 19800000250109301 1 98330 1 1 9800000 
430Q9 86O1 1894 3010340 1 1 9046QQ9660 1 2QQ31 Ol 1 9001 19125Q1 1 9flnn«5f57 
1201092000014600894013002601 197018264901 10602601 19804 1973301 
190000001 101092000014601 126014QQ15Q1 189000001 101 19800005430n 
98601 1894401094010933201 197000002601 18901092260000001 1974900 
0000482601 151 0046926005170000049O0 5O6ononnnnnonoooonnnnnnnnn 



* * oooooooooooooooooz 

1 600 5360 1 256 1 5p 1 75300002 1 60 1 1 1 20QQ99320 1 1 900000Q49004580430 1 
27601 19Q490096604301290O1 1 £9482201 18901 168250132501 1763301 17 
600000150 1 168 0000 1 6p 1 34760(K)0 1 101347Q0Q01 21 01 17601 17643Q139 



3101 380038584903826 
3101 440038584903826 



4011 68490 1 3500 160144101 83^50 1 4280 1 1 69 1 1 1 44 1 *00230 1 1 75000 
00260 1 477Q 1 44 U gQ 1 477Q0QQ4260 1 1 66000002201 1 6600093320 1 1 ft?ono 



3101 500038584903826 
310 1 560038584903826 
31016200 38584 903826 
31 1 680038584903826 



002301 17601 1652601 186019372201 1 86000963201 1 780000023Q 1 1 860 1 1 
86 3 30 00 B£ QGO-Q 1.5 008 1 oooo J2 1 o oo9on 1 1 Q6P6n l tflftnn n 9n23ni i ftAni i 
652 1 0009500095 1 20 1 3470000 1 470 1 6220 1 200260 1 1 8600094440 1 6940 1 3 
253201 185 00 Q00230 1 1 8 6Q 1 1 97 1 1 p 1 1 89QQ QNQ250 1 093Q0Q99260 1 0920 1 1 



3101 740038584903826 
3 1 0,1 800038584903826 



3101860038584903826 
3101 9200385£4903826 



892601 19800092460177401300490096604700882014004901 1260000000 
000000O0OOOOOOOOO00OOOOQOOOnJ9n6QgRnflnJ74n7TP^ftJ6niPflQAAR.]A 



3101 980038584903826 
31 02Q40038934_9_03826 



8249N490 J380 1 3M759 J29099M 1 64 J2 1 267L674 Jl 4332L266 J081 47K922 JO 
2597K630 J POOOOQO OO 1 6Q Q 5360 1 970 1 50 1 753000034Q0 1 2240430 1 QQOO 1 1 
9049009660430201001 169490096603201 169000001201 16800GN12301 19 
* 7011762101 1890116 S4!9Q17 



I 
I 



2 



31021 14038584903826 4800000000003703759005002503769004003 1 00048o3758360383800500 
31 02174 03858490 3826 2 1 02 196Q384 1 1 6QQQ64jQ39502600P79n3fl42 1 1 n3B4?nnnKnPfi0n0690.1R4? 
31 02234038584903826 2600074038474302326038482 1 022880006436038380050031 0000 1 03838 
31 022 94038584903826 1 1 02288000Q047022700Q90049000 1 203603643005 00 1 403645000004603 
3102354038584903826 1 740 1 2002303648000692 1 00099000693300099000002 1 00098036452 1 00 
310241 40 385849 03826 09800 063330364600 0003 2 0009 5000 002603 160000994 4Q2490036482602 
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I # PROGRAM DESCRIPTION 



This program solves linear programming problems with 
output of detailed results* That is, given coeffi- 
cients a£4, cost coefficients 03, and requirements 
bi, determine xj such that 

£ a ij*3 - bi with xj^O 
and 

E. clxi * maximum 
J ° 

Computations are performed by the Dual Algorithm until 
a feasible solution is obtained, control is then given 
to the simplex Algorithm for optimization, cost changes 
and requirement changes can be made after loading 
original matrix or after solving original matrix. 



4. Machine requirements 

!• Basic 1620 with card input and output. 

2. Any size storage can be used. The larger the 
storage, the larger the problem that can be 
solved. The size of the problem is restricted 
by the following relationship: 

(m +2) (n + 3) ^ memory - 3920 



where; m is the number of restrictions 

n is the number of non-basis independent 

variables 
memory is 20,000, ij.0,000, or 60,000 

3. source language: actual machine language 

subroutines: floating point subroutines supplied 
with program. Ho other subroutines needed. 

5. Input media: card reader and console typewriter. 

6. Output media: card punch and console typewriter. 

7. Running time: The precise time required per 
iteration depends on the size and density of 
the matrix. As an approximation, a problem 
with 30. equations and lj.0 non-basis variables 
requires about 20 seconds per iteration. 

8. Decimal accuracy: All computations are performed 
in 2-and-8 floating point form* Matrix input 
can be either fixed point or floating point. 
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B« Error stops 

1* All stops may bo recognised by displaying 
IR-1 in MAR. MA.R will contain the address 
of the halt instruction plus 12, 

felt ins tr« Reason 

01126 Floating point overflot* 

01288 Floating point attempted division 

by zero 

00012 Loader - ready for next code 

022*58 cost changer - Change entered for 

a non-existent identification 

number 

00012 cost changer - ready for next code 

03110 RHS changer - change entered for 

a non-existent identification 

number 

00012 RHS changer - ready for next code 

025 3k Dual algorithm - inconsistent matrix 

0265k simplex algorithm - unbounded solution 

O2838 Matrix punch - ready for next code 

2. Error procedure 

a* Cost changer ( 021*58 )• Erroneous identification 
number is the last one typed* 

If data are being entered through the type- 
writer, depress RESET and INSERT. 
Type lj.9 02126. 
Depress RELEASE and START. 
Corrected or new change data may now be 
entered normally. 

If data are being entered from cards, the 
erroneous datum may be ignored by following 
the same steps as above (for entry through 
typewriter). 

If it is desired to enter the change correctly 

through the typewriter, the following steps 

must be taken: 

Set console switch two on. 

Depress RESET and INSERT. 

Type U9 02126. 

Depress RELEASE and START* 

Bater correct datum through typewriter 

(complete with record mark). 

Depress RELEASE and START. 

!4hen typewriter is again selected, depress 

RESET. 

Set console switch two off. 
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Depress INSERT, RELEASE, and START. 
Processing will now continue normally, 
reading data from cards* 

b. RHS Changer (03110): Erroneous identification 
number is the last one typed. 

If data are being entered through the type- 
writer, depress RESET and INSERT. 
Type k9 02126. 
Depress RELEASE and START* 
Corrected or new change data may now be 
entered normally. 

If data are being entered from cards, the 
erroneous data may be ignored by following 
the same steps as above (for entry through 
typewriter). 

If it is desired to enter the change correctly 

through the typewriter, the following steps 

must be taken: 

Set console switch three on. 

Depress RESET and INSERT* 

Type 1^9 02126. 

Depress RELEASE and START. 

Enter correct identification number (ij. digits 

plus record mark) through typewriter. 

Depress RELEASE and START* 

Enter requirement change (including record 

mark) through typewriter* 

Depress RELEASE and START* 

When typewriter is again selected, depress 

RESET. 

Set console switch three off. 
Depress INSERT, RELEASE, and START. 
Processing will now continue normally, reading 
data from cards. 

c. Dual algorithm (02£&): BJo obtain solution 
existing at the time of discovery of inconsis- 
tency: 

Run out cards from reader. 

Replace cards in reader hopper starting with 

Pinal Basis Output program. 

Depress READER START. 

Depress RESET and INSERT* 

Type ij.9 0002l|.* 

Depress RELEASE and START* 

d* simplex Algorithm (02652$.): To obtain solution 
existing at the time of discovery of unbounded 
solution: 



o 



o 



Depress RESET and INSERT. 

Type JU-9 0002!^, 

Depress RELEASE and START. 

e. Overflow-underflow: Overflow causes a 

stop at 01126 or 01288, Underflow causes 
the result to be set to zero, and process- 
ing continues* Either condition indicates 
*i$her machine malfunction or unreasonable 
data* 

f • Any other error: unrecoverable • 
Restart problem. 



II, THE PROGRAM 



A, Operation instructions 

1, Set typewriter margins at 12 and 92; tabs at 
2^, 3o, k$t 60, and 69, 

2, It is important to note that no more than one 
instruction at a time may be inserted at any 
time after the initiation of the Iioader routine, 

3, To load initial data: 

a. Ready the LOADER deck in the card reader 
(including matrix input data). Depress 
READER START. 

b. Depress RESET and INSERT. 

c. Type 16 00010 00000. 

d. Depress RELEASE and START. 

!£his initiates memory clear routine. 

e. After at least one second, depress INSTANT 
STOP, RESET, and INSERT. 

f • Type 36 03826 00500 
k$ 03826 

g. Depress RELEASE and START. 

h. Program deck loads, data are loaded, and a 
halt at location 00012 is executed. Opera- 
tion may now proceed to the COST CHANGER, 
RHS CHANGER, or SOLUTION programs. 
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!}!• To effect cost changes: 

a. Ready the COST CHANGER deck in the card 
reader (including data change cards if 
changes will be made from cards). Depress 
READER START. 

b. set console switch two off if cost changes 
will be entered from card reader; switch 
two on if data will be entered through the 
typewriter. 

c. Depress START. 

d. If data are entered through the typewriter, 
RELEASE and START must be depressed after 
each entry. (Dontt forget the record mark 
to terminate each entry.) 

e. After all data are entered a halt at location 
00012 is executed. Operation may now proceed 
to the RHS CHANGER or SOLUTION programs. 

5. To effect requirement changes: 

a. Ready the RHS CHANGER deck in the card 
reader (including data change cards if 
changes will be made from cards). Depress 
READER START. 

b. set console switch three off if requirement 
changes will be entered from card reader; 
switch three on if data will be entered 
through the typewriter, 

c. Depress START. 

d. If data are entered through the typewriter, 
RELEASE and START must be depressed after 
the entry of each identification and each 
change. (Don't forget the record mark 

to terminate each entry. ) 

e. After all daba are entered a hs.lt at location 
00012 is executed. Operation may now proceed 
to the COST CHANGER or SOLUTION programs. 

6. To solve and obtain output: 

a. Ready SOLVER deck in the card reader. 
Depress READER START. 

b. Set console switch one on to type iterations, 
off to omit iteration type out. 
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c. set console switch four on to obtain puncii- 
out of final matrix, off~To suppress punch- 
ing • 



d. Depress START. 



e. After all output is complete, a halt at 

location 02838 will be executed • Operation 
may now proceed to the LOADER, COST CHANGER, 
or RHS CHANGER programs. 



B. Data Preparation 

1. The first card to be entered is a case identification 
card* A five digit alphabetic Identification 

is punched In columns 1 through $• columns 6 
through 80 are blank. 

2. l*he second card specifies the size of the problem 
being entered* It is punched as follows: 

Information card columns 

(zero with a minus sign) 1 
number of equations (3 digits) 2 - If 

(zero) $ 
(zero with a minus sign) 6 
number of non-basis variables (3 dig.) 7-9 
(zero) 10 
input code 11 

a. 1 for row- column, fixed point 
input 

b. for a complete floating point 
matrix 



Columns 12 through 80 are blank. 

3. UJie matrix ean be entered in one of two forms: 
floating point (the entire matrix is entered* 
column by column, exactly as it is to be stored* 
including zero entries) or fixed point (entries 
can be in any order, accompanied by their res- 
pective row-column designations, and zero ele- 
ments need not be entered). 

Por floating point entry, eight elements are 
punched on each card, with a minus sign over 
the first position of each ten position field, 
and a minus sign over the last position of 
negative fields (units position). There are ten 
positions per field, with the first two positions 
specifying placement of the decimal point (50 
designates a decimal point to the left of the 
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first non-zero digit; 51 designates one non- 
zero digit to the left of the decimal; 14.9 desig- 
nates that the first non-zero digit is in the 
hundredths position, etc.) and the last eight 
digits (with no leading zeroes) giving the 
element itself. A zero element Is specified 
by ten zeroes, with a minus sign over the first. 

The elements must be punched in column order, 
starting with the first basis cost* See the 
matrix layout included with this writeup to 
know how the elements are stored in the matrix. 

!?♦ Fixed point matrix elements are entered four 
to a card. 

Two types of cards are necessary for matrix 
elements : 

a. row- column cards designate where the elements 
will be stored in the matrix. Row-column 
designations are -punched in columns 1-6, 
21-26, il-ij.6, and^ 61-66 in the form RRRCCC 
(a minus sign accompanies the first digit ^ 
of each row and column designation). The 
column designation specifies column number 
considering only non-basis vectors* The 
requirement variables have only a row desig- 
nation, and 000 should be entered for their 
column. 

b # Matrix element cards contain fixed-point 
entries including optional leading sign, 
up to eight numeric digits, and mandatory 
decimal point. Each entry must be terminated 
by a record mark (0-2-8 punch). The four 
elements are punched beginning in columns 
1, 21, I4.I, and 61 for each card. 

Elements may be entered in any order, but corres- 
ponding row-column designations and matrix 
entries must be in the same positions of suc- 
cessive cards (row-column card preceding). 

Loading of zero elements Is optional 



If any card contains fewer than four elements, 
the first blank field (column 21, Ijl, or 6l) 
of the row-column card must be punched with a 
record mark (0-2-8 punch). 
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Item Description Row-column 

a^j matrix element of the ith row -^c^ 
and jth non-basis vector xxxxxx 

h± element in the ith row of 

the requirement vector xxx<?00 

6m Fixed point cost entries are also four per card, 
but the row-column designation is on the same 
card as the cost. 

Entries consist of a six digit row-column desig- 
nation followed by a ten digit ID/cost, where 
the ID is a four digit numerical designation 
to identify the cost variable, and cost is a 
six digit fixed point entry with the decimal 
after the third digit* 

Item Description Row- column 

ci cost per unit of basis <— — * v 

variable for row i xxxTToTxxxxxxxxxx 

J X"D <^s +- 

c* cost per unit of 3 th _ ^^T* — v/ — * ^ 

non-basis vector OOTxxxxxxxxxxxxx 

Negative cost is signified by a minus sign over 
the units position of the cost field. 

If any card contains fewer than four elements , 
the first blank field (column 21, [p., or 61) 
of that card must be punched with a record mark 
(0-2-8 punch). 

7. i*he_last data entry for fixed point mode toast 

be 000 to terminate the loading operation* #his 
can be punched in columns 1-3 of a nevr card, 
or in 21-23, ip.-ij.3> or 61-63 or the last row- 
column or cost entry card. 

8. Entry of cost changes: 

Data for the COST CHANGER routine can be entered 
either from cards or from the typewriter. If 
entered from cards, one entry is punched in 
each card, in columns 1-10. The entry consists 
of the four digit ID and six digit cost figure, 
with a minus sign over the first digit of the * 
cost (column 5>) and over the last digit of negative 
costs (column 10). Following the last cost 
change card must be a card punched with four zeroes 
(columns to terminate the operation. 
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If entered from the typewriter, the format is 
the same as from cards: four digit ID and six 
digit cost figure with flag over the first digit 
of all cost figures, and over the units position 
of negative costs. The last entry must be four 
SSeroes to terminate the operation. 

9» Entry of requirement changes: 

Uiere are two entries for each requirement change 
the identification and the change in the value 
of that requirement element. Note that the 
new requirement element is not entered, but 
rather the difference between the old and new 
requirement elements. If changes are entered 
from cards, the identification and change are 
entered on separate cards, one per card. 

The identification is entered first, and con- 
sists of the ID associated with the cost element 
for that row. Thus, to change bi ■ b3, use the 
ID associated with C3. = 03. Ho flags or record 
marks are needed in entering the 3D either 
from cards or typewriter. It is punched in 
columns 1-lf of card entries, with $-80 blank. 

After entering the ID, its associated requirement 
change is entered (by typewriter or on a separate 
card). The entry starts in column 1, if a 
card is being used, and consists of optional 
leading sign, up to eight numeric digits, and 
mandatory decimal point. It must be terminated 
by a record mark, whether it is entered from 
cards or from the typewriter. 

The last entry must be four zeroes to terminate 
the operation. 

Interpretation of Results 

1. Cost changes or requirement changes are typed 

if the input is from cards. If input is through 
the typewriter, the input operation itself 
produces a log of the changes made. 

2. Dual Algorithm and simplex Algorithm: 

Results of iterations are typed only if console 
switch one is on. If it is desired to monitor 
the course of the solution onlh occasionally, 
this may be done by setting switch one on and 
off periodically. The following points should 
be borne in mind: 
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a # Pag© headings for iteration output are 
typed only if switch one is on at the be- 
ginning of the solution* 

b. The iteration count is reset when the dual 
algorithm is finished and control passes 
to the simplex algorithm. 

The information which may be typed includes 
the iteration count, the current value of the 
profit function (floating point), the last 
variable to leave the basis, and the variable 
which entered the basis. 

Final Basis Output; 

The final value of the profit function is given 
in floating point. 

For each final basis variable, the following 
information is supplied (in fixed point): 

a. Identification/cost coefficient. 

b. Activity. 

c. The limits of the cost coefficient over 
which the current solution is optimal. 

d. The variables which limit the range of the 
cost coefficient and will enter the basis 
if a limit is exceeded. 

Final Non-Basis Output: 

For each non-basis variable, the following 
information is given (in fixed point): 

a. Identification/cost coefficient. 

b. shadow price - the penalty to the total 
system if a unit of this variable is forced 
into the final solution. 

c. The limits of activity over which the shadow 
price applies. 

d. The variables which limit the range of 
applicability of the shadox* price* The 
upper limiting variable is the one which 
would leave the basis if the associated 
non-bg.sis variable were forced into the 
solution. 
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A punchout of the complete final matrix is op- 
tional under control of console switch four." 
The cards punched also contain all necessary 
supplementary information and are in a format 
which is suitable for direct re -loading by 
the LOADER routine. 

' • Functions of Individual Sub-programs 

1m Uie data load routine reads and stores parameters 
for the problem. It then proceeds to load and 
store the matrix in the proper format. In the 
case of fixed point input, storage locations 
are computed from the row/column designations, 
and the input elements are converted to floating 
point notation as they are read and stored. 

2. The cost change routine reads ID/cost elements, 
searches the matrix for a corresponding iden- 
tification number, and stores the new cost 
coefficient in the proper location. 

3» Hh.e requirement change routine reads variable 
identification numbers and changes for the 
associated requirement element. The change is 
converted to floating point notations, the 
matrix is searched for the identification num- 
ber, and the requirement vector is updated. 

ij.. The shadow price calculator converts all cost 
coefficients to floating point notation and 
evaluates Zj-cj for each non-basis variable. 

$m The dual algorithm computes a feasible solution 
for the problem. Its operation may be monitored 
allowing the typing of the value of the functional 
and the identifications of the variables that 
enter into and are removed from the basis. It 
should be noted in this regard that the objectives 
of feasibility and optima lity may not be compati- 
ble at this point, and therefore the functional 
may actually decrease from one iteration to the 
next while the dual algorithm is in control. 

During the reduction of the matrix of coefficients 
a test is made to determine whether the resulting 
quantity has a magnitude less than an amount 
called "essential zero." Any element falling 
into this category is set to ssero in order to 
avoid computation with numbers which are actually 
not significant. "Essential zero" has been set 
to 10"°. Its value (the corresponding floating 
point exponent) is located in memory addresses 
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311j.ll. and 31^5 and may be changed if desired* if 
this change is made, it is important to note 
that the high order digit must carry a flag. 

6« The simplex algorithm starts with any feasible 
solution and computes the optimal feasible solu- 
tion* Its operation may be monitored by allow- 
ing iteration typeout as in the case of the dual 
algorithm. 

l!he discussion of essential zero as applied to 
the dual algorithm in 5. above is equally valid 
for the simplex algorithm. In this case the test 
value is located in memory addresses 32% and 
3265. 

m the case of a tie in the quotient used as the 
criterion for choosing the variable to leave the 
basis on any iteration, the choice is resolved 
by using the largest denominator as a secondary 
criterion. 

7» The basis output routine searches out and com- 
putes the information described in section c. 
above. Floating point quantities (with the ex- 
ception of the functional) are converted to 
fixed point before being typed. 

8. The non-basis output routine searches out and 
computes the information described in section 
C. above. Floating point quantities are converted 
to fixed point before being typed. 

9m The matrix punch routine resets the permanent 

instructions in low memory to allow re-initiation 
of the input or change phases. Problem parameters 
and the final matrix are punched into paper tape 
if called for by the setting of console switch 
four. 

Use of storage 

Memory locations from 00012 through 03919 are used 
by the program. Memory from 03920 up to the re- 
quired location is used for the matrix and its mani- 
pulation. Locations 00000 through 00011 are avail- 
able for the insertion of one instruction at a time 
through the typewriter. A halt instruction in posijrf 
tion 00012 ensures that only this one instruction 
will be executed. 

00012-0001^ load routine 
OOOL0- 00058 case identification 
00060-0006ij. s address 



-13- 



00065-00069 M and 2 

O0070-0007l|- H 

00075-00079 M 

00080-00099 product area 

00100-00i}.00 arithmetic tables 

OOlj.02-02061^ floating point subroutines 

variable programs 

03758-03917 input area (record mark in 03918) 

03920- matrix 



Program Decks 
1. LOADER 



3 pre-load cards. 
Arithmetic tables, 
boating point routines. 
Data load routine. 
Matrix data to be entered. 



2. COST CHAKSER 



2 pre-load cards. 
Arithmetic tables. 
Floating point routines. 
Cost changer routine. 

Cost change data (if entered from cards). 



3. RMS CHANGER 



2 pre-load cards. 
Arithmetic tables. 
Floating point routines. 
Requirement change routine. 

Requirement change data (if entered from cards). 
SOLUTION 



2 pre-load cards. 

Shadow price calculation routine. 

Page headings. 

2 pre-load cards. 

Dual algorithm routine. 

2 pre-load cards. 

Simplex algorithm. 

2 pre-load cards. 

Final basis output routine. 

2 pre-load cards. 

Fix and print routine. 

Page headings. 

2 pre-load cards. 

Non-basis output routine. 

Page headings. 

2 Pre-load cards. 

Matrix punch routine. 
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a* Storage of matrix in memory 

The matrix is stored in column sequence and uses 
memory from 3920 upward to the extent required by 
the problem size, oSie following diagram indicates 
the manner in which matrix elements are stored: 



3920 






S4l0(2m43 ) 3D/C j S+10[(nfl ) (m*2 




S 


Functional 


S410(2mf l±) S+10 £(n4l ) (m42 )J 


Working 
Column 

39204 
lO(mfl) 


ID/Ci 
S*llQm 


SU0(m*3) 
S410(2m+2) 


S410 (2m4$ ) 341o[W 1 ) (m+2 )«.] 

Increasing J 

a lj5 S410|~(541)(m+2)4iJ 

Inc ceasing 
i 

S+10(3m+i}.) S*10{jnf 2) J-2)-2] 



S is computed for each problem and is equal to 3920 
4 10(m43). All addresses are field addresses* 
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III. SAMPIE PROBLEM 
A. Matrix 

Requirement 

Basis variables Non-basis variables vector 



9901 


9992 


0001 


0002 


9902 




1 




1 


2 




6 




I 


1 


-1 


-1 







-10 


1 


1 








V&riable 
ID 



Cost 



B* Matrix in storage 





0000000000 


S001001000 S002001000 9902000000 




0000000000 


Sooooooooo Sooooooooo ffOOOOOOOOO 


5*901000000 
999201000(5" 


5i6ooooooo 
5U4.0000007 


5110000000 5L20000000 0000000000 

5110000000 ?iiooooooo £110000000 



2X3 
£0020^00301 
001001 

002703 

hlm£ „ 

TTOlt»OT2901000000 
S0lff039902000000 

Coo 



Input 

!• Fixed point 



2X3 

00020606300 



S01$02 
2.* 

Soiffoo 

^SBcBBsisfr 01000$ 



2« Floating point 



tf02ff01 
ff02j00 

T&iJoiSooiooiooo 



S02S02 

SoTOo277oo200iooo 



NiinoonoooNiiooooonoooo2ooioooooooooooooNi2oooooonNnooooonnR9n2oooooooooooonndft 
Simon. _^r> Nllonnnnnnnnoonnoooonooonoooooonoonoon n oono o oon o^o^nn. v>oo r H l Hvn n nnn i i f ^HinV 



© 



o 
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Output 

1. Pinal matrix punched 



2X3 

00020000300 



^i3ooooooQ N120000000R992 

2. Typewriter output 



CASE 

ITER NO 

501 
502 
001 

FUNCT IONAL 

VAR/COST 

5902000000 
0001001000 



2X3 

FUNCTIONAL 

0000000 
513OO0OO00 
5160000000 
J160000000 



VAR OUT 



VAR IN 



599201 0005 0002001 000 
9901 000000 5902000000 
0002001 000 0001 001 000 



ACTIVITY LIM VAR 



10.0000 
6.0000 



5002 
5002 



VAR/COST SHAD PRICE LIM VAR 



5002001000 
5992010005 

5901000000 



1.0000 
10.0000 

1.0000 



LOWER LIM 


LIM VAR 


UPPER LIM 


.3333- 


5992 


10.0000 


.5000 


0000 


9999.9000 


LOWER LIM 


LIM VAR 


UPPER LIM 


9999.9000- 


5ooi 


3.0000 


10.0000- 


5ooo 


9999.9000 


9999.9000- 


5ooi 


6.0000 



COST CHANGEB^- 



f READ IN 
ARITH. 



IV, METHOD OP LOADING PROGRAM 

The "LOAD" key cannot be used to load programs in this 
routine since the first 80 positions' of memory are not 
free. Therefore cards are read into locations 03838 
through 03917> with a record mark in 03918 being utilized 
to transfer up to six instructions to their place in 
memory. The first two instructions on each card con- 
tain the "transfer record" instruction and a branch to 
location 03826, which contains the read instruction 
36 03838 00500. The first two cards of each routine 
perform the function of setting th» read instruction 
in 03826 and the record mark in 03918. The first card 
is read into 03826, and control is transferred there, 
using instructions 36 03826 00^00 I4.9 03826. This 
must be done by the operator ail the beginning of the 
"LOADER" program, it is done automatically for all 
follow! ng routines. 




READ NEW 
COST FROM 
TYPEWRITER 



READ NEW 
COST PROM 
CARP S, i 



5r 



LACE NEH 

COST IN 
MATRIX 



RHS CHANGER 



SHADOW PRICE 



© 
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MANUALLY 
CLEAR 

MEMORY 



GENERAL PLOW 
Page 1 



READ IN 

ARITH* 

TABLES 



jfffiAD" 
FLOATING PT 
SUBROUTINES 



fREAD IN 

LOADER 

(PROGRAM 



READ IN 
MATRIX DATA 



FLOAT & 
STORE IN 
MATRIX 



RHS CHANGER 



SOLUTION 



feAD IN j 
ISHADOW PRIC 
[PROGRAM I 



■ FLOAT 
PRICES 
P MATRIX] 



"11 



[READ IN 
PAGE 

Ljheadi 



"1 




(READ IN \ 
DUAL ALGO. ; 
1 PROGRAM i 



COMPUTE } 
FEASIBLE \ 
SOLUTION : 





READ NEW 
REQ. FRO* 
TYPEWRITER 



(READ NEW 
REQ. FROM 
CARDS 



t> r <] 



PLACE )l£\i 
REQ. IN 
MATRIX 
/ 



SHADOW PRICE 



\ 

CoIt CHANGER 
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GENERAL PLOW 
Page 2 




[COMPUTE 

[optimum 
Feasible 




[READ IN 
jfON-BASIS 0, 
PRO GRAM 

f READ IN 
PAGE 

READINGS 



PAGE 
HEADINGS 



TYPE 

NON-BASI0 
OUTPUT 



/READ IN 
MATRIX PCH. 
[PROGRAM 



PUNCH _ 
MATRIX Sf 
SW l± ON 



COST CHANGER 



LOADER 



RHS CHANGER 
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V 22SS 

INITIALISE 
STORE 
INSTR. 



5 




DATA LOADER 
Page 1 



r 33?6 
| READ 
| ROW/COL 

L 




r 



2270 
READ FL. 
INPUT 
CARD 

STORE 

IN 
MATRIX 

f 229ty. 
MODIFY 
PREV. 
INSTR* 




£ — 



(IE 
317fr 



COMPUTE 
MATRIX SIZE 
10(N+2)(M42V 



CARD 

O- 

>w=0 v ;^-(no) 



COMPUTE 4 STORS 
STORAGE LOCATION 

s jio[(j»i)-(i»2hi ] 




© 



o © o 
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_Jt-2210 

SET 
;J0BK AREA 
TO ZERO 
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211b. 



SET i=0 
SET x-ra 



SHADOW PRICE CALC, 
Pag© 1 



2150 



2l5 




i = i+1 

MOVE 0^ 
TO 

WORK AREA 




,2198 
'Ci«0?V-/ HO) 





SET 


EXP 


■ 53 


i ^ 


= l 


$2290 




Formalize] 

PILL OUT 
TRAILING- 
ZEROES 





SHADOW PRICE CALC. 
Page 2 



X - 3C-1 




Pg. 1 



MODIFY" LOOP 
TO FLOAT & 
SBSATB Cj's 



SET j-0 
SET x«n 



7 a go 



3 8 3U 



MOVE Ga 
TO J 
*ORK AREA 




£2^0 



SET 




EXP - 


53 


k = 


i 



¥2290 



SAVE 
SIGH 




k»k*l 
fcXP*EXP- 



U3- 



239k 



SAVE k 



X 



', FORMALIZE 



PILL OUT 
TRAILING 
ZEROES 
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© 



SHADOW PRICE CALC. 
Page 3 



SET 
OPPOSITE 
SIGN 



STORE FLOATED 

3 1 in shadow 
Price row 



^2tt90 




SET 
FUNCTION ih 
TO ZERO 



[SET J»0 



Y27S2 



SET i=l 
SET y=m 



-^27^ 
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SHADOW PRICE GALC. 
Page k. 




TYPE 
CASE 
NUMBER 



READ 
PAGE | 




V pit 

READ 
NEXT 
CO PE 

I 

I 

DUAL 
ALGORITHM 
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2134 



SET 


= 


SET x - 


m 


SET i - 






17= 




141 ! 



?!7k 



^ES) <5^n? — ^o) 

_.$2198 
iJET buiin 




/ NO) c X « OtV-YYEl 



NEXT 




CODE 





SET W 

SET y»n 
SET j=0 



SIMPLEX 
ALGORITHM 



iSET ROW 4 

Address 1 



J = J+l 



DUAL ALGORITHM 
Page 1 



2330 





DUAL ALGORITHM 
Page k 



SET j=0 
SET y=n 



O 

no62 




^5 a U 

,3 _=r— 

y = y-i 




PIPE: ITERATION, 
FUNCTIONAL, 
ED/cOSffc 
ID/COST r 



T 
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S IMPLEX ALGORITHM 
Page 1 



( 2126) 



\ 




SET dmin = 


SET y 


s n 


SET j 


. 1 



J 2210 t 

SET dj^sdj j 

i t 



J = i+1 

y = y-i 





SIMPLEX ALGOfilTHM 
Page 2 



SIMPLEX ALGORITHM 
Page 3 



12810 



MOVE 



SET a^k 
- 



i = 141 
x = x-1 



f HO 



882 




:<0? 



28% 



SET 
ark « 1 



SET j=0 
SET y=n 



\}?%? 



a rJ* a rj/ a rk 




SET i«€ 
SET J»0 

SET 




SIMPLEX ALGORITHM 
Page ij. 



1122 



SET j«0 
SET y«n 



VnfiP 
| a 13' a ir ( a P3 )(a l jjj 




3706 



j 


= 3+1 


y 


= y-i 




SWITCH 
ID r & IDi; 



STEP 
fcTERATIOl* 

COUNT 



TYPE: ITERATION, 
FUNCTIONAL j 
ID/COSTfc 
ID/COST r 




Pg. 1 
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2562 



READ 
FIX 
UB ROUTHfe 



FINAL BASIS OUTPUT 
Page 1 



TYPE 

:otctiona|l 



HEAD & 

TYPE PG< 
1EEADINGS 



Jim 



SET i=l 
SET x=m 



X29Q6 



3ET UPPER » -Oo 
5ET LOWER = «*> 
SET J ■ 1 
SET y - n 




'2906 






pAVE 











f 321k 


3 * 

y * 


7-1 



3190 



' SET 








(SAVE 






o 



o 



FINAL BASIS OUTPUT 
Page 2 



mo 



TYPE TD/C± 



3fXPE ACTIVITY] 



tSTE LOWER ITH 



ty%06 

SET C=Ci| 



E 

& type! 



GjALC 
LOWER LIMIT 
C- (LOWER) 



t 



|fype itpp^r' nj 



t 



CALC. & TYPI 
UPPER LIMIT 
• C- (UPPER) 



Jj63P 



1 = i-l 
X = x-1 



s3*78 



fYEa 



rx = o?; 



[NO) 



2906 



Pg. 1 



L2t 



HEAD 
NEXT 
CODE 



V 

NON-BASIS 
OUTPUT 



O 



TURN ON 
ALT FLAG 



X 



FIX SUBROUTINE 
Page 1 



SET EDIT 
WORD IN 
CUTPUT ARIA 




FORM RT 
JUSTIPIE1 ' 
FIXED PT 










SET 






FIAG 






i 


— 












ADD 


c 1 



SET WORD 
TO 9 f s 




FIX SUBROUTINE 
Page 2 



^EsV~<^i 



&ETBLANK 
IK OUTPUT 

digit 3 : 

%>288 




i = 141 



. _ V2276 , 

MOVE DIGIT 
1 TO OUTPUT 
pIGIT J 



-J - J-H 



/SJl-^ -sc - o? >. /no",_ 




2598 



NON-BASIS OUTPUT 
Page 1 



TYPE 
PAGE 
HEADINGS 

I 



SET j=l 
SET y»n 



SET UPPER- o° 
SET LOWER= -°° 
SET 1=1 
SET x=m 




2802 




"UPPER?>-^ES) 



SET 
UPPER 
= R 



{sAVE IDj] 





$3110 



i = i-1 




NO 



0.9k 
0> 



t^206 



HPE U>/C j 



TYPE 
SHADOW 
PRICE 



© 



-ilk- 

^ NON-BASIS OUTPUT 

Pago 2 




LOWER ID 



TYPE 
LOVER 
LIMIT 

_Jf 

TXPE 
UPPER I ) 



TYPE 
TIPPER 
LIMIT 




IttTRIX 



MATRIX PUNCH 
Page 1 




HALT 



READ 
NEXT 

CODE 



LOADER 



I RHS 
^ CHANGER 



2610 



PUNCH 
CASE 



_2&70 



PUNCH 
«,N, & 
CODE 



COMPUTE 
MATRIX 
LIMITS 



T 2778 



PUNCH 
CARD OB 
MATRIX 




COST 
CHANGER 



PRE-LOAD" 
CARDS. 



ARITHMETIC 
TABLES 



360383800500 
1 503918000024903826 
31 00012038874903826 
31 001 00038584 90 3826 
31 001 60038584 90 3826 
31 002200385^84903826 
31 00280038564903826 
31 00340038584903826 
3100^02038584903826 
3100462038584903826 

31 00 522038584 903826 
31 00582038584 90 3826 
310064 203858490 3826 
31 00702038584903826 
3 1 00762038584903826 
31 00822038584903826 
31 00882038584 90 3826 
31 0094203858490 3826 
310 1 002038584903826 
310 1 06203858490 3826 

3101 122038584903826 
3101 182039004903826 
31 01200038584903826 
31 01260038584903826 
31 01 32 0038584903826 
310 1 380038584 90 3826 
31 144003858490 3826 
310 1 500038584903826 
310 1 560038584903826 
310 1 620038584903826 
31 1 680038584903826 
31 01 740038584903826 
310 180 0038584 903826 
3 1 1 860038584903826 
310 1 920038504903826 
3101 98003858490 3826 
31 02040038934903826 



* 4800000000003603826005004903826 

00000000000 1 0203040 00204060800 030609021 00408021 61 00500 1 51 02 
0060218142007041 1 28200806 1 422300908 1 726300000000005060708090 

1214161815181 1242720242822363520353045403632484455324946536 
04646^462754453627 180 1 ^34^6789 i 234367890234567890 J34567890JK 
4567890 JKL 567890 JKL^67890JKLMN7890JKLMNO89OJi<LMNUP90JKLfMNOPQ 
160053600 53849004340160053600558260 1 19700O0o260l 1 12004451601 
12400^001 lOl 1240000 1460049401400490 1 1280046926O1 176000002601 
18701827490000002201 18601 1764900570 02601 18601 176150084700001 
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16030970 393026034400007922039300 1 826460321801200260308502221 
1603121000002603416000742603 14003121 1 203 1 4000008 1 60046903Q97 
2601 1970 0000490 1 938039301 600469031331 60044500000490040200099 
1400000000^3470359801 30021 030850006921 0314 00006921 031 2100069 
12034 160 00 J 4603062013001 102984000 JO 1 103097000J01 1 03025000 JO 
1 203440000 J 460297801 300260337302221 26H3361 O35601 203361 000 JO 
260336803361 2203373000692603380033732600023O0 000260000000000 
26000000 023 1 1 03631 00001 47021 14 001 0034nnnonno 1 0238o362900 1 00 

340000000 1 08260357703373260354 1 03361 26035050006421 0350500069 
26039 1 700000380390800 1 00340 000000 1 0826039 1 700000380390800 1 00 
34 000Q000 1 08260391 700000380390800 1 0049021 1 4n26036 1 603 1 2 1 260 
* ■ . * 0000182649031 56000 



25036680 374726000 2301826260354800074260219700064210219700069 
21021970006924000230000047022340 1 100260223302 197260002300 000 
2 102,1 970 0069 120354 8000 JO 4702 1860 120 04400024O002326000230 3522 
1603481 00000 260363200 079260238002233260243302233260242 100064 
21 02421 00 69490254202400 000 01 826470251 801 1 00 1 60046902433260 1 
1 9700000490 1 20 000 000240009900023460 25 1 80 1 1 004 6025860 1 2 002602 
49302380260 374600000260002300 099260 34810 2421 1 203632000 JO 4 602 
63001200 1 102380000 JO 1 102421 000 JO 1 1 02433000 J04902 374026026090 
23802403746000004 70247001 300490251801403481 0000047026660 1 200 
48N1 10000000 260280 40 24 93260 28280 2233260363200079260278502493 
1 60281 60 3930260290002493260282 1 02233 1 60046902785260 1 1 9702665 
490120000000260392000099260000001826260393000000260000001 826 
1 1 028 1 60 00 JO 1 102821 000 JO 1 1 02828000 JO 1 203632O00 J046028 1 00 1 300 
26000000266526029770348126029840 3481 260354800074 1 60046902977 
2601 1 970 3920490 1 9380 0OO0260 0000000992 1 02977r»006921 02984Q0069 
1 203548000J0 460294201 30016031 0403930260314500064210314500069 
1 60321 703930260363200079220^9300 1 V^w46o3338o 1 20026032050348 1 
160324 10 0000 2603548000742603*2600324 1 1 20326000 008 1 600469032 1 7 
260 1 1 970 0000 490 1 93803930 1 600^6903253 1 &00445o000049004020O099 
1400000000M347037060 1 3002 10320500069210 32600006921 03241 00069 
1 2035480 0OJO4603 1820 13001 103 104000 JO 1 1 0321 7000 JO 1 1031 45000 JO 
1 203632000 J046030980 1 300260346902233 1 2o3469oOO J0260347603469 
220348 1 00 69 26 O 34 880 34 81 2 6Q OQ23 OQ 00 026o000000000260000000023 
1 10366700001 47021 26Q0 100 340000000 102380366500 100 34Q0000001 08 
26036850 3481 260364903469260361 30006421 03,61 300069260374600 000 
380373700 1 00 340000000 1 08260 374600 00038o3737o0 1 00 3400000O01 08 
26037460 0000380373700 1 004902 1 26026037240324 1 26000000 1 8264903 
* * 2780000000000Z 
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PRE-LOAD 
CARDS 



PINAL 
BASIS 
OUTPUT 
ROUTINE 



PAGE 

HEADINGS 



360383800 5QQ 
-1 50391 80000Z4903826 
'31 02562038584903826 
31 02622038584903826 
31 0268203d584903826 
31 02742038584903826 
3 10 2802038584? 03826 
31 02862038584903826 
31 02922038584903826 
31 02982038584903826 
31 03042038584903826 
31 03102038584903826 
3103 162038584903826 
31 03222038584903826 
31 03282038584903826 
31 033420 3858490 3826 
3 1 03402038584903826 
31 034620 38584 90 3826 
31 03522038584903826 
31 03582038584903826 
31 03642038584903826 
31 03702038904902562 
360383800500 
1 5039180000Z4903826 
3 1 01 94 8038584 903826 
31 02008038584903826 
31 02068038584903826 
31 02128038584903826 
3102188038584903826 
31 02248038584903826 

3102308038584903826 
3102368038584903826 
31 02428038584903826 
31 02488038584903826 
3102548039044903826 
2503919004004902582 
FUNCTI ONALZ 

VAR/COST Z 
ACTIVITY Z 
LIM VAR Z 
LOWER LIM Z 
LIM VAR Z 
UPPER LIM Z 



3603826005004903826034001 00001 02260 
75900500390 375900 1 00260391 70000034 
0000010234 000 0001 02 1 603533000 06370 
5330000 147027140 1200 3400000 001 02260 
369026531 1 03369000 J026030 13033692 1 
1 000J02602965030 1 31 1 02965000K0 260 
740025931600 72500000160071500000160 
1 09029651 60308500000 16031 5700000260 
000 03706460 321 40 1 200 16004690309 7260 
1 66000 240009900 750470321 4ol 1 00260 
2 1 4024000990074046032 1 40 1 3002600740 
006921030850006921 030970006921 03109 
00691 203533000 JO 47030 380 1 200260 3918 
0100 320 1 955000 00 2701 948000 00260*91 8 
1 026034290 3321 26007600000 3 *0 0755 

0000 320 3521 00000 490 34980 330 3521 0000 
330 1 9550 00 00 3200 74 0000002701 9480074 
380391 50 01 OO 3300 75000000 3301 955 00 00 

1 10332 10 00J0 1 1 03369000 JO 1 203545000 J 

* 000 



2653000642 1 02653000693703 
001 0001 08L803908001 0034 00 
3759005003903759001 001203 
3321000641 103321000J02603 
3013000 792603001030131 1 03 
3545000792600 750026662600 
30440 O0OO2603097029652603 
321 3030 132603533000742200 
1 197000004901200000004403 
07500 00^92600 73000000490 3 
000992600 7200000021 030440 
0006921 03-1 570 00 6921 0321 30 
000^0 34 000 0^00 1 02380 390 9 O 
007143400000001 063S0391 50 
0000044034860 07 603300 7600 
0260076900760 1 50077000000 
02603918007243400000001 08 
027O1948007501 1 02965O00 Jo 
04 7O2906O 1 20049 02400000 
0070 707070037070707000070 



320 19560 00003 1038980370744020080 1947330 194700 000 16O391700OKO 
26007600 3706140 1 939000N4470206401 1 00 1 600757R99994902 1 480 1 20 1 
939000M74702 1360 1300 1602 1350 19402 102 1350 19391 50 193P0000021 00 
7600194744024420 1955160221900004 1602287007531602282038991602 
243007531602250038993201 9570000044022760 19574302264007531603 
7590000049022880330195700000250389900 7531 i 0228700001 1 1 022430 

000 1 1 10 2282000 02 1 1 0225000002 1 20 22 1 90 000 1 4 70 222 00 1 20 044024 1 60 
1956330 195600000330 195700000 16022 1900 00549 023 1 20 340000000 108 
3903899001 0042140 391 7000K04 7024 780 120032007600000021 00760007 
7047025220 1300160391 7000004902 1480 1 6039 1 7OOOK04902 1 48707O707 
* * 00370 707070000 
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PAGE 
HEADINGS 



PRE-LOAD 
CARDS 

MATRIX 

PUNCH 

ROUTINE 



360383800500 

1 5039 1 80000Z4903826 

3102562038584903826 34000 1 000 1 02L400 1 0000 1 02 1 60325300006370375900500390375900 1 00 
3102622038584903826 1 2032530000 1 47025980 1 200340000000 1 0226032 1 7000642 1 032 1 700069 
31 02682038584903826 2 1 032 1 7000692603265032 1 71 2032 1 7000 J0260286 1 032 1 7 1 1 02861 OOOKo 
3102742038584903826 26028970286 1 220289700069260290902897220290900069260322400074 
31 0280 2038584903826 26007500 2585260074002574 1 6007 1 500000 1 60072500000 1 60294000000 
3102862038584903826 26029930286 1 26030050286 1 1 60298 100000160 305300 000260310902909 
3102922038584903826 26032530007922000000343846031 1 00 1 200 1 6004690^993260 1 1 9700000 
31 02982038584903826 490 1 20 00 000044030 6200000240 00990 074 04703 1 1 no 1 1 00260074000099 
3103042038584903826 26007200000049031 1 0024 000990075046031 1 00 1 3 00^600750000992600 
3103102038584903826 730000001 102940000J01 102981000J01 102993000J01 103005000J01 103 
3103162038584903826 053000 JO 1 1 03 1 09000 JO 1 203253000 J0470 29340 1 20026039 1 8000003400 
31 0322203858490 3826 0000 1 02380390900 1 00320 1 95500 00270 1 9480000026039 1 80 07 1 4 3400 
31 03282038584903826 00000 1 0838039 1 500 1 002 70 1 9480074026039 1 800724340000000 1 083803 
3103342038584903826 9 1 500 1 00270 1 94800 7502 1 0286 1 0692 1 32 1 70 00 692 1 0326500069 1 203 
31 0340 20388 1 4902562 * 224 00 J04 70 280 20 1 20 049000240000000000 

VAR/COST Z 
SHAD PRICE Z 
L I M VAR Z 
LOWER LIM Z 
LIM VAR Z 
UPPER LIM Z 
360383800500 
15039180000Z4903826 

31 02562038584903826 4702826004003 1 035630340 1 3 1 0356200048 1 50357300000390356300400 
31 02622038584903826 3 1 035630 348 1 2603567000792603572000 74 1 5035 7 300 003803563004 00 
31 02682038584903826 2602784 064 1 1 0278400 00 1 230 0069000733200 095O0 0O02 1 0009900069 
31 02742038584903826 2 1 00098000782 1 000990278426028 1 3000993800000004 00 1 1 02784000Q0 
31 02802038584903826 1 4 027840 000047027780 1 3004 1 0000000 J0048888888888849000240000Z 

3JL9340 1 038584903826 — - : : — — 

3t,p346 1 038584903826 : =■ 

3103521038784903826 * — — 

250356 1004002503723004 04902562 
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TECHNOLOGY 



BBWIl 



4% 



4 



4 



4 



4 



7" 
bv 



N 

4 



7 

bv 



V 



A K 



4 



7 




4 



4 



bv 



A 



4 



V 

K 



4 



4 



17 
bv 



7 




4 



bv 



V 



4 



4 



7 
bv 



4- 

"7 



7 



4 



bv 4 



7 ^ 



b. 



V 



\ A 

7 ^ 



4 



bv A 



V N 



7 
bv 



^"7 



F 



4 



^ A 

7 ^ 



bv 



7 



4 



4 



\ 



bv 



V 



7 



A N 



4 



4 



^ 4 



F N 



bv 



4 



7 ^ 



V 



4 



7 



bv 4 



F 



F 



4 



4 



bv 



A A. 

-7 



A. 



7 
A 

7 



"7 



4 



bv 4 



7 N 



1 

7 
bv 



17 N 



7 



4 



bv 4 



b* 



4 



7 



4 



17 \f 



bv 4 



4 



^1 



V 
bv 



4 



bv 4 



7 ^ 



4 



bv 4 



17 

bv 



V 



A N 



bv 4 



7 x 



4 



4 



7 ^ 



|7 



4 



bv 4 



^1 



v 

bv 



7 



A N 



bv 4 



7 N 



4 



7 ^ 



2_ 



a. 
V 



V 
7 



4 



bv 4 



7 N 



7 

bv 



"7 



bv 



7 



7 



bv 4 



bv 4 



7 \ 



4 



7 



7 
bv 



4 



bv 4 



7 N 



7 
bv 



4 



7 



bv 2! 



7 ^ 



17 
bv 



4 



4 



17 ^ 



7 



K 4 



7 ^ 



C7 



4 



7 ^ 



4 



K 4 



V 
bx 



t7 



4-b. 



b. 4 



TP 7 



4 



7 ^ 



-7 ^ _ 



A. 



7 



t7 
A 

7 




^1 

4 



b. 4 



b. 



bv 



^ A 
7 ^ 



V 



A N 



b. 4 



\TT7 



4 



b. 4 



A 



7 ^ 



4 



V 



7 



t» 4 



K 4 



P 7 N 



4 



17 



b, 4 



V ^ 



^1 
4 



bv 



w 

A N 



K 4 



7 



4 



.bv 



A A. 



7 



PUROUE UNIVERSITY 



CALUMET CA* 



V N 



4 



7 



7 ^ 



V 



7 ^ 



4 



P 7 



bx 4 



4 



bv 



4 



bv 4 



V 



^ 4 



4 



bv 



A 



7 ^ 



bv 4 



7 



4 



F 

bv 



bv 4 



F N 



7 



iv A 
7 \ 



F 

bv 



V 

bv 



bv 



F 



A A. 

-7 



4 



bv 



7 



bv 4 



7 N 



4 



bv 



iv A 
7 ^ 



4 



V 

bv 



bv 4 



7 \ 



4 



F 

bv 



b. 4 



7 N 



4 



F 

bv 



bv 4 



V N 



4 



t7 

bv 



7 



bv 4 



7 ^ 



4 



bv 



N 

4 



bv 4 



F 

bv 



F 



7 



h 4 



F N 



4 



bv 



4 



4 



7 ^ 



7 

bv 



F 



A. 

-7 N" 



7 



4 



4-bv 



bv 4 



F N 



bv 4 



f 



n~F 



4-bv 



\~7 



bv 4 



F N 



4 



bv 4 



7 ^ 



N 

4 



7 

bv 



bv 4 



7 N 



4 



7 

bv 



bv 4 



7 N 



4 



7 

bv 



7 



bv 4 



7 'q 



4 



bv 



bv 4 



4 



7 



'q 



bv 



7 



7 



bv 4 



7 N 



4 



bv 



iv A 
7 \ 



4 



bv 



7 



7 

bv 



N~7 



bv 



7 



■7 ^ _ 



4-bv 



"7 



v 



7 

bv 



bv 



F 



F 



bv 4 



F ^ 



4-bv 



bv 



bv 4 



7 N 



n~F 



4-bv 



F 



bv 



4-bv 



iv A 



V N 



^"7 



7 



7 



7 ^ 



4 



bv 



7 \ 



4 



7 



7 

bv 



bv 



7 



7 



bv 4 



4 



bv 



bv 



7 



'q 



7 

bv 



bv 



7 



4 



F 

bv 



4 



F 



F 

bv 



N 
A 



bv 



V 

7 



C7) 

M 

55 

n 

g 

o 



td 

g 

3 



I ^ 

o 
o o 

a* 

•i O 

S (o 



1620 Correction 
10.1.006 



Feb. 1,1965 
0101 



' "Linear Programming Code for the Card 1620 vlth Punched card 

Option for Final Output" 



An error in the subject program was reported to the 1620 Users Group 
and subsequently to us. The enclosed changes should eliminate the 
error. 

The Right Hand Side Changer Program (cards numbered 200 - 260) should 
be modified to eliminate a stray flag bit left in memory position 03513. 
The following changes should be made. 

1. Card #229 of the RHS changer program. Change cols. 1-7 from 
4800000 to 4903484. 

2. Card #259 of the RHS changer program. Change from col. 33 on to 

contain: _ «. 

330351300000480000000000*000001510345203508300259 

3. Card #260 of the RHS changer program. Change cols. 1 - 13 to 
490212600000* and 63 - 80 to 

010350803520000260 

4. Add Card #261 to the RHS changer program. Card cols, and contents 
as follows 

Col 1 - 69 - all_zeros 
Cols 70 - 80 - 02114000261 

5. Appendix A, Section D3 should read: 

"RHS changer itself loads and a stop at 03507 is 

executed." 



A complete corrected RHS Changer deck has been distributed since 
Feb. 2,196"5 
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File No. 10.1 



LINEAR PROGRAMMING CODE 

FOR THE CARD 1620 
WITH PUNCHED CARD OPTION 
FOR FINAL OUTPUT 



Lou Davis and Art Nickel 
IBM 

401 Grand Avenue 
Oakland, California 
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l£20 Correction leb. 1,19*5 

10.1.006 0101 

"Linear Programing Code for the Card 1)620 with Punched card 
Option for final Output 11 



An error in the subject program was reported to the 1620 Users Group 
and subsequently to us. The enclosed changes should eliminate the 
error. 

The Right Hand Side Changer Program (cards numbered 200 - 260) should 
be modified to eliminate a stray flag bit left in memory position 03513. 
The following changes should be made. 

1. Card #229 of the RHS changer program. Change cols. 1-7 from 
4800000 to 4903484. 

2. Card #259 of the RHS changer program. Change from col. 33 on to 

contain: 

330351300000480000000000*00000310345273508000259 



3. Card #260 of the RHS changer program. Change cols. 1 - 13 to 

490212600000* and 63 - 80 to 
010350803520000260 

4. Add Card #261 to the RHS changer program. Card cols, and contents 

as follows 

Col 1 - 69 - all.zeros 
Cols 70 - 80 - 02114000261 

5. Appendix A, Section D3 should read: 

"RHS changer Itself loads and a stop at 03507 is 
executed." 



A complete corrected RHS Changer deck has been distributed slnot 
Feb.2,19^5 



1620 10. 1. 006 



DECK KEY 

1. Data Loader - 000-053 

2. Cost Changer - 100-107 

3. BBS Changer - 200- 260 

4. Solution Deck - 300-916 (not continuous) 



A. 


Floating Point Subroutine 


300-328 


B. 


Shadow Price 


400-423 


C. 


Dual 


500-526 


D. 


Simplex 


600-628 


E. 


Final Basis 


700-720 


F. 


Fix and Print 


750-768 


G. 


Non Basis 


800-821 


H. 


Matrix Punch 


900-916 



5. The deck mentioned in Section I. on Page 10 



© 



3. 



PROGRAM ABSTRACT 



Title: Linear Programming Code for the Card 1620 with Punched Card 
Option for Final Output 



Subject Calssification : 10.1.006 



Authors : Lou Davis and Art Nickel 
IBM 

401 Grand Avenue 
Oakland 10, California 



Direct In quiries to : Lou Davis or Art Nickel 
IBM 

401 Grand Avenue 
Oakland 10, California 
TEmplebar 4-7070 



Purpose : Solution of linear programming problems with output of 

detailed results. Given coefficients a, j, cost coefficients 
Cj , and requirements b t , determine xj such that 

^^i # j x j = bj with xj *2?0 

and 

^c-x. = maximum. 

i J J 



Mathematical Method: Computations are performed by the Dual Algorithm 
until a feasible solution is obtained. Control is 
then given to the Simplex Algorithm for optimization. 
Many, many things go on before this stage is reached 
and after. It is quite important to read the instruc- 
tions for order of program input (Appendix A) and data 
input carefully. 

a. Accuracy: All computations are performed in 2- 
and -8 floating point arithmetic. 

b. Derivation -Reference: Some (Nichols') notation 
and techniques were derived from the writeup of 
the "Linear Programming Code for the Augmented 
650" by O.R. Perry. Reference is also made to 
C. R. Nichols' writeup for the 1620 paper tape 
input/output version. 



© 



Restrictions, Range : 

a. 



b. 
c. 



Storage Requirements : 
Equipment Specifications : 



Source Language: 



Program Execution Time : 



Check Out Status: 



Requires a 1622 Card Read-Punch Unit. This 
program was rewritten for a 20K machine. Certain 
changes in the program deck are necessary to 
enable it to run on a 40K or 60K machine. These 
changes are indicated in Appendix E. The size 
of the problem which can be handled is restricted 
by the following relationship: 

(m+2) (n+3) ^ memory - 3920 
10 

where m is the number of restrictions, 

n is the number of non-basis independent 
variables, and memory is 20K, 40K, or 60K. 

Data must be prepared in the format specified in 
Appendix B . 

Output may be either on the typewriter or on cards. 
The optional final matrix punchout is on cards. 
(See Addendum No. 1 to program writeup). 



Any size memory - see Restrictions. 



Basic 1620 with Card input and output. 



The original tape program and the modification 
for card input and output have been coded in 
1620 SPS. 



The precise time required per iteration depends 
on the size and density of the matrix. As an 
approximation, a problem with 30 equations and 
40 non-basis variables requires about 20 seconds 
per iteration. 



This program has been widely and successfully 
used in the field for some time now on many prob- 
lems of different sizes and descriptions. 
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COMMENTS 



This program and its documentation were written by an IBM 
employee. It was developed for a specific purpose and submitted 
for general distribution to interested parties in the hope that it 
might prove helpful to other members of the data processing 
community. The program and its documentation are essentially 
in the author's original form. IBM serves as the distribution 
agency in supplying this program. Questions concerning the 
use of the program should be directed to the author's attention. 
For certain problems it has been found that the value for "essential 
zero" located in memory positions 3144 and 3145 may need to be 
increased to some value greater than the 10-8 value to which it is 
set in the program. 
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I. General Information 

II . Appendix A 
HI, Appendix B 

IV. Appendix C 

V. Addendum No. 1 

VI. Appendix D 

VII. Appendix E 

VIII. Addendum No. 1 
DC. Addendum s 



- Operator Information 

- Input Format 

- Output 

- Card Output of Final Basis 
and Non-Basis 

- Sample Problem 

- Storage layout and listings 

- Punched Final Output 

- Overflow, Negative Data and 
Conversion of "Punch Option" 
Deck for 40 K or 60 K. 



1620 UNEAR PROGRAMMING CODE FOR CARD INPUT/OUTPUT 



Nichols, Nickel and Davis 
Addendums to 1620 linear Programming Code - File Number 10. 1. 006 



To avoid an overflow condition and certain loading inefficiency when using 
the Data Loader program, make the following change in the program deck 
and listings for the Data Loader: 

Object Program 

Card No. Location From To 

32 01946 4902022 4902090 

Symbolic listing (Data Loader) 

Label Location From To 

TSTRC + 12 01946 B Fill . B Signst 

Ihe subject "overflow" causes no apparent error but only slows down 
the execution of the program. If the above change is made, data cards 
must not contain record marks. 

A second item that it would be well to note is the following: 
When using negative (minus) data values, the minus sign (11 punch) 
i must be in the first and only first column of the ten (10) column field 
allowed for the number. That is, column 9, 29, 49 and/or 69. 
If this is not followed, the value will be treated as positive. 

To convert the "final output - punch option" deck for use on a 40K or 60K 
1620, the following cards and columns should be changed in the manner 
prescribed in the program writeup: 



t 



o 



Card Number Columns 

700 3 and 51 

701 3, 15, and 39 

702 27 and 39 
712 47 

718 43 

800 39 and 51 

814 23 

912 15 

These columns contain ones (l r s) as the first digit of 19xxx addresses 
and should be changed to threes (3's) or fives (5 ! s) for the 40K or 60K 
machine. 
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1620 Linear Programming Code for Card Input/Output 
Nichols, Nickel and Davis 
Addendum No. 1 to Program write-up 

I. Reference: Item 3. "Restrictions", parte. 

The Final Basis and Non-Basis Output may be punched on cards 
rather than typed by replacing cards numbered 700 through 916 of 
the solution deck with a different card deck bearing the same first 
and last number (700 and SI 6) and approximately the same numbers 
in between. (These cards are the last section in the deck. ) 

II. Reference: Appendix C - Output 

Items 2 and 3 

2.1 and 3.1 

Final Basis and Non-Basis Output using "Punch Deck" for cards 
700 - 916. The final value of the profit function is typed in 
floating point form as before. 

The same information is punched as is given with the typed output. 
The card format is as follows: 

Card Columns 

1 - 2 Code 11 - Final Basis 
22 - Non- Basis 

3-10 Zeros 

11-20 Identification/Cost Coefficient 

21-24 Zeros 

25 - 32 Activity - Final Basis 

or Shadow Price - Non- Bask 

33 Record Mark 0-2-8 punch 

34 Zero 

35 - 38 Limiting variable 

10 



39 - 40 Zeros 

41 - 48 Lower Limit on Cost - Final Basis 
or Lower Limit on Activity - Non-Basis 

49 Record Mark 0-2-8 punch 

50 Zero 

51 - 54 Limiting variable 
55 - 56 Zeros 

57 - 64 Upper Limit on Cost - Final Basis 

or Upper Limit on Activity - Non- Basis 

65 Record Mark 0-2-8 punch 

66 - 80 Zeros 



Notes: 1. The card output is numeric so there are no decimal points 
or minus signs. 

2 . Negative numbers are denoted by a flag (l 1 punch) over the 
low order digit i.e. , card columns 48 or 64. 

3. The Activity, Shadow Price and Limits are 8 digit numbers, 
and the decimal point is assumed in the middle — that is, 4 
positions from the right or left, so (xxxx.xxxx). 

4. Card columns shown as containing zeros may contain "most 
anything" on the first card punched, but all columns of any 
interest will be correct. 



1. Identification : 

a. 1620 Linear Programing Code for card input /out put. 

b. this program is an amended version of C. P. Nichols* 1620 
Linear Programing Code for paper tape input/output. Moat of 
his original program and much of his writeup have been retained 
unchanged for this edition. This version edited by Art Nickel 
and Lou Davis. 

c. IBM Branch Office, Oakland, California. 

2. Pur pose : Solution of linear programming problems with output of 
detailed results. Given coefficients ai,^ cost coefficients cj, and re- 
quirements bj, determine Xj such that 

£ a i,j x J = b i vithxj^O 
and 

^ CjXj^tiaximuaw 

j 

3. Restrictions : 

a. Requires a 1622 Card Read -Punch Unit. This program was rewritten 
for a 2 OK machine. Certain changes in the program deck are necessary 
to enable it to run on a 4 OK or 6 OK machine. These changes are in- 
dicated in Appendix K. The size of the problem which can be handled 
is restricted by the following relationahip; 

(m + 2)(n + 3)£ —"frf 3960 
where m is the number of restrictions, 

n ia the number of non-basis independent variables, and 
memory is 20K, UoK, or 60K. 

b. Data must be prepared in the format specified in Appendix B. 

c. All output ia on the typewriter except an optional final matrix 
punch out which la on cards. 
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k. Method : Computations are performed by the Dual Algorithm until a feasible 
solution is obtained. Control is then given to the Simplex Algorithm for 
optimization. Many, many things go on before this stage is reached and after. 
It is quite important to read the instructions for order cf program input 
(Appendix A) and data input carefully. 

a. Accuracy: All computations are performed in 2-and-8 floating point 
arithmetic . 

b. Derivation-Reference: Some (Nichols') notation and techniques were 
derived from the writeup of the "Linear Programming Code for the 
Augmented 650" by 0. R. Perry. Reference is also made to C. R. 
Nichols' write up for the lb20 paper tape input/output version. 

5. Miscellaneous : 

a. Special Features: 

(1) Cost changes in terms of the original data can be made without 
the necessity of reloading or resolving the original matrix. 

(2) Requirements can also be changed — as in (l) above. 

(3) The overall program consists of several relatively independent 
subprograms which can be deleted or altered in many respects without 
interference with other subprograms. Sequential loading of the 
subprograms is automatic, with intermediate programmed stops at 
points where it might be desireable to manipulate order of input 
decks. 
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APPENDIX A 
OPERATOR INFORMATION 

A. Preliminary : 

1. Set typewriter margins at 15 and 95; tabs at 27, 39, 51, 63, and 72. 

B. Load matrix : 

1. Clear Memory « turn off all console switches. 

2. Place DATA LOADER deck followed by data deck in 1622 Card Reader. 

3. Make sure computer la in manual mode. 
h. Depress LOAD on card-reader. 

5. DATA LOADER itself loads and a halt at OOU13 is executed. (All 
addresses given for halts will be the address displayed in the 
Memory Address Register. ) 

6. Depress START on console. 

7. If all goes mil the data will load and a normal halt at 02265 
will be executed. 

8. Operation may now proceed to the COST CHANGER, RHS CHANGER, or 
SOLUTION. 



C. Cost Changes : 

1. Place COST CHANGER deck (followed by coat changes if input ia by 
cards) in card-reader hopper. 

2. Press START on console and on card-reader. 

3. COST CHANGER itaelf loads and a stop at 02125 i« executed. 

fc. a. If coat change* are on cards make sure Sense Switch 2 is off 
and depress START on console. 

b. If coat changes are to be entered from the typewriter, turn 

Sense Switch 2 on. (IELBASE mad START must be depressed after 
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aach entry—don't forget last entry of 000C^# ) 
If all goes well data will load and program will branch to the sub- 
routine to load the next program in the card hopper. If there Is no 

program in the hopper, the computer will hang up on a READER NO FEED 
* 

with a MAR address of 19851 * In this case put the appropriate deck 
in the hopper and push START on the reader. 

It is possible for all not to go well because a change has been en- 
tered for a non-existent identification number. In this case the 
computer will halt at 02** 57 . 
a. Error Procedure: 

(1) If data are being entered through the typewriter, depress 
RESET and INSERT, type U902126, depress RELEASE and START. 
Corrected or new change data may now be entered normally. 

(2) If data are being entered through the card reader, NON PROCESS 
RUNOUT cards In hopper, find offending card, correct it, put 
cards back in hopper at that point (i.e., the corrected card 
followed by cards yet to be processed), press START on card 
reader, depress RESET and INSERT, type 4902126, depress 
RELEASE and START. 

Alternatively, after card(s) has been corrected it is possible 
to reload the entire deck, Including the COST CHANGER deck, 
and go through the procedures in (2) except type I+9198U0. 
Operation may now proceed to RHS CHANGER, or SOLUTION (or COST CHANGER 
again ) . 

Right Hand 8ide (Requirement) Changes ; 

Place RHS CHANGER deck (followed by RHS change data deck if Input is 
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by cards) in card reader hopper. 

Press START on console and on card reader. 

RHS CHANGER itself loads and a stop at 02125 is executed. 

a. If RHS changes are on cards make sure Sense Switch 3 is off 
and depress START on console. 

b. If RHS changes are to be entered from the typewriter, turn 
Sense Switch 3 on, and press START on console. (RELEASE and 
START must be depressed after each entry— don* t forget last 
entry of 0000$) . 

If all goes well the data will load and the program will stop at 
03U75. 

Just like the COST CHANGER it is possible for strange happenings 
because a change has been entered for a non-existent identification 
number. In this case the computer will halt at 03201 . 
a* Error Procedure: 

(1) If data are being entered through the typewriter, depress 
RESET and INSERT, type 1*902126, depress RELEASE and START. 
Corrected or new change data may now be entered normally. 

(2) If data are being entered through the card reader, NON PROCESS 
RUNOUT cards in hopper, find offending card, correct it, put 
cards back in hopper at that point (i.e., the corrected card 
followed by cards yet to be processed), press START on card 
reader, depress RESET and INSERT, type 1*902126, depress 
RELEASE and START. 

No alternate, unl aaa y e w. t a k e- w% o a rd e a lr ea dy p reeeaaed i 

(Oee ab ote e, 5., a., (£ ! )>. ■ 
Operation may now proceed to COST CHANGER or SOLUTION (or RHS 
CHANGER again). 
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E, Solve the problem ; 

1. Place SOLUTION deck in card reader hopper. 

2. Set Sense Switch 1 on to type iterations, off to omit iteration 
typeout . 

3. Set Sense Switch h on to obtain punchout of final matrix, off to 
suppress punching. 

U. Depress START on console (and if needed on card reader). 
5. Error Stops: 

a. 025U5 Dual Algorithm-Inconsistent Matrix 

b. 02665 Simplex Algorithm-Unbounded Solution 

It may be possible to obtain the solution existing at the time 

of either of these stops by NON PROCESS RUNOUT the card reader, 

put back in the hopper that portion of the SOLUTION deck be- 

# 

ginning with card number ''00, RESET and INSERT, type U91?BU0, 
RELEASE and START on the console and START on the card reader. 
b. After all output is complete, a halt at location 028U5 will be 

executed. Operation may now proceed to DATA LOADER, COST CHANGER, 
or RHS CHANGER. 

F. Other Stops : 

01137 Floating Point Overflow 

01299 floating Point Attempted Division by Zero 

While overflow causes a halt at either of the two above addressee, 

underflow causes the result to be set to zero, and processing continues. 
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APPENDIX B 
INPUT FORMAT 



Prologue: There are three basic types of data input: 
Input for DATA LOADER (A below) 
Input for COST CHANGER (B below) 
Input for RHS CHANGER (C below) 



A. Input for DATA LOADER : 

1. First Data Input Card - Header Card. The first card of the data 
deck must be a card containing the case identification. Only the 
first five columns (columns 1-5 ) are read by the program. Explicitly, 
columns 1-5 of the first data card may contain any alphanumeric 
characters to identify the problem being run. 

2. Second Data Input Card - Parameter Card. The second data card muat 
contain certain Information needed by the program. Namely, 

Columns Information Example 

1-3 Number of restrictions. 006 
U Blank 

5-7 Number of non-baaia variables 008 

8 Blank 

9 Input type code (one digit) as follows: 

1 for row-column, fixed point input 
for a complete (include zeros) 
floating point matrix 

3. She data, following the above two cards, can be in one of two formats: 
a. Row-Column, Fixed Point Input: Up to four elements with row-column 

designations may be entered on each card for this type of input. 
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In general, the format is as follows: 
i j element 
XXXbXXXbXXXXXXXXXX 
where b stands for blank column. 
There are four fields per card for data: 







Columns 
for i 


Columns Columns for 
for J element 


Field 


1 


1-3 


5-7 


9-18 


Field 


2 


21-23 


25-27 


29-38 


Field 


3 


U1-U3 


45-47 


49-58 


Field 


4 


61-63 


65-67 


69-78 



All other columns are not used by the program and may be used for 
anything else. 

ANY field which has the middle digit of the three digit i value left 
blank will be ignored and the other columns of that field may be used 
for anything else. Hence, data cards do not have to have a complete 
set of four elements, nor is it required that the left hand fields 
be used first, nor is it required that consecutive fields be used. 
— > ffatrix elements may contain up to eight numeric digits with optional 
leading sign and mandatory decimal point. Loading of zero elements 
is optional i 

IMPORTANT: The signal to the DATA LOADER to quit loading is a three 

digit value of i-000. This must be in the last data card in a field 
the 

toA*lght (if other values on the card) of the laat element. The best 
thing is to have a card with 000 punched In columns 1-3 and then have 
this card as the last card in the data deck. 
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(Fixed Point Input continued) 
Item Description 



Row -Column 



Matrix element in the ith ^ 
row and of the jth non-basis XXX XXX 
vector. 



Remarks -Examples 

010 oo4 36.075 

003 004 -4.0 

008 017 -.6 

008 017 -0.6 



Element in the ith row of 
the requirement vector 

Jr 000 

Cost per unit of the basis 

variable for row i 
j- 001= OOJ where 
J is a flagged one. 



XXX 000 



XXX OOJ 



008 000 6.025 



Expressed in fixed point 
with the following format: 
ID COST 

xxx ooj IcxaJfiocxxx' 



(1 



Decimal point in coat assumed 
here, negative cost denoted 
by flag (11 punch) over units 
position of cost. 

012 OOJ 0003001500 Here 
cost is 1.5 

012 OOJ OOI6OO5999 Here 
coat is -5.999 



Cost per unit of the Jth 

non -basis vector 

i — 001 ■= OOJ where 
J is a flagged one. 



Remarks same aa above, 

OOJ 024 0034123456 Ha 
coat is 123- ^56 



Summary of DATA LOAD input for row-column, fixed point input. 




/XXXXX 
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b. Input for Floating Point Matrix: In the floating point mode, the 

entire matrix (including zero elements) is loaded in column order. 

(See storage layout in Appendix E). Seven ten digit floating point 
matrix 

numbers must be punched in/column order in each card in columns 1-70. 
Elements of a new matrix column follow immediately the last element 
of the old matrix column (i.e., a new matrix column does not 
automatically start a new card). A record mark (0,2,8) punch must 
be placed in column 71 of each card. If the last card does not 
require a full 7 elements, the record mark must instead follow 
directly behind the last digit of the last element. ALSO, the 
i ast card must have a record mark in column 72 (or else you will not 
terminal* loading). Columns 73-80 you have for your own pleasure. 
The floating point representation is the same as used in lb20 SPS. 
Note that the optional matrix punch out of the SOLUTION deck pre- 
pareB a dat*deck in this format. 



8-2 X 
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APPENDIX C 
OUTPUT 

1. Results of iterations are typed only if Sense Switch 1 Is on. If it 
is desired to monitor the course of the solution only occasionally, this may be 
done by setting Sense Switch 1 on and off periodically. The following points 
should be borne in mind: 

a. Page headings for iteration output are typed only if Sense Switch 1 
is on at the beginning of the solution. 

b. The iteration count is reset when the dual algorithm is finished 
and control passes to the simplex algorithm. 

The information which may be typed includes the iteration count, the 
current value of the profit function (floating point), the last variable to leave 
the basis, and the variable which entered the basis. 
?.. Final Basis Output: 

The final value of the profit function is given in floating point. 

For each final basis variable, the following information is supplied 
(in fixed point): 

a. Identification/cost coefficient. 

b. Activity. 

c. The limits of the cost coefficient over which the current solution 
is optimal. 

d. The variables which limit the range of the coat coefficient and will 
enter the basis if a limit is exceeded. 

3. Final Non-Basis Output: 

For each final non-basis variable, the following infarction is given 
(in fixed point) : 

L 
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a. Identification/cost coefficient. 

b. Shadow price - the penalty to the total system if a unit of this 
variable is forced into the final solution. 

c. The limits of activity over which the shadow price applies. 

d. The variables which limit the range of applicability of the shadow 
price. The uper limiting variable is the one which would leave the 
basis if the aaaociated non-basis variable were forced into the 
solution* 

$. A punchout of the complete final matrix is optional under control of 
Sense Switch h. The deck punched out contains the Header and the Parameter 
cards and is in a format which is suitable for direct reloading by the 
DATA. LOADER routine. 
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B. Input for COST CHANGER ; 

1. Cards: One element only per card. For each cost coefficient to be 
changed, a ten digit ID/new cost element is punched in columns 1-10 
in the following format. 

ID New Cost 

Decimal point assumed here. ' 
Negative cost is signified by a flag over the units positions. 

2. Typewriter: As above— ten digits typed in without any spacing 
between digits, but immediately followed by a record mark . 

3. IMPORTANT: Termination of loading is recognized by a field of four 
zeros, 0000. These are punched in columns 1-k if entry is 

or typed in as OOOOtif entry is by typewriter. 



C. Input for RHS CHANGER : 

For each requirement element to be changed, a four digit identification 
and the change amount are entered. The identification entered is that for 
the original basis variable for the row whose requirement is to be changed. 
The identification format is XXXX. 

The second entry for each change is the amount bjr which the original 
requirement is to be clanged . Its format is XXXXXXXXXX. 

1. Cards: One element only per card. 

Columns 1-4 Columns 5-8 Columns 9-18 Example 

ID Blank Change Amou nt 

AJtXX bbbb WOOOOCXXXX' 0012 -0.75 

The change amount may be any number (fixed point) up to eight digits 
with mandatory decimal point and optional leading sign. 



if 
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but a record mark must immediately follow each entry. 
Typewriter. Sane as above/ RELEASE and START will have to be 

depressed between entry of the ID and the Change Amount. 
IMPORTANT: Termination of loading is recognized by a field of four 
zeros ^ 0000. These are punched in columns 1-k if entry is by card 
or typed in as 0000*1 f entry is by typewriter. 
Cost changes and/or requirement changes say be made in unlimited 
numbers on any single run. If both types of changes are to be made, 
the order In which the two changers are operated is of no consequence. 
It should be noted that changes may be initiated under the procedures 
given in Appendix A either immediately after loading s matrix or at 
the conclusion of the output phase. If the output phase is inter- 
rupted before its completion, however, the computer is not conditioned 
to accept change b. 



Sample &<>8leaa 

f t f IX ^ - Tnin 




tlx + n +y fc 



= 4 
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TEST1 
005 004 1 

001 001 1. 001 002 3. 001 003 -1. 002 001 2,0 

002 002 1.0 002 00* -1.0 

003 001 -1. 003 002 3. 

004 001 1.0 

004 002 1. 

005 001 1. 
005 002 -1. 

001 000 6. . 



003 000 ,. ^ ^jr^>r>> 

004 000 6. A* ^ ■ » -^\(4i* 1 



005 000 2. 

001 OOJ 000599999R 002 OOj 000699999R 003 OOJ 0007000000 004 OOJ 0008000000 
005 OOJ 0009000000 

OOJ 001 0001001000 OOJ 002 0002002005 OOJ 003 0003000000 

00 OOJ 004 0004000000 

000 



-rw pe^nvz/L output 5a*apl£ c**>e^ 

CASE TEST1 



ITER NO FUNCTIONAL VAR OUT VAR IN 

001 5420039980 0005999999 0002002000 

002 $144000000* 0006999999 OOOIOOIOOO 



3«, I , 



FUNCTIONAL 


5 144000000 










VAR/COST 


ACTIVITY 


LIM VAR 


LOWER LIM 


LIM VAR 


UPPER LIM 


0002002000 


1.6000 


0004 


3.0000- 


0003 


.5000- 


0001 00 1005 


1.2000 


0003 


4.0001- 


0004 


.6667- 


0007000000 


6.4000 


0003 


.4285- 


0004 


.1666 


0008000000 


3.2000 


0004 


.5000- 


0006 


2499.^975 


0009000000 


2.4000 


0004 


.2500- 


0003 


1.0000 


VAR/COST 


SHAO PRICE 


LIM VAR 


LOWER LIM 


LIM VAR 


UPPER LIM 


0OO6999999 


999.7990 


0009 


3.0000- 


0001 


2.0000 


0005999999 


999.3989 


0007 


4.5714- 


0002 


4.0000 


0003000000 


.6000 


0O02 


4.0000- 


0007 


4.57* 


0004000000 


.2000 


0001 


2.0000- 


0009 


3.0000 



D 5 



0009 1. 
0OO4 -1»5 
0000 



OdOlOOlOOO 
0002002000 

oooo 



APPENDIX E 
STORAGE LAYOUT AND LISTINGS 

I. STORAGE USEUr 

Memory locations froir 00048 through 3747 are used by th- program. Locations 19840-19963 are used for e program 
lo?d routine. Locition^ IP9f>4, ISS65-00044 are used for a crd in;ut area for data, heMings, and program input (it 
is # Viwover, possible to usi 00000-0^0047 for inserting). 'Aemory from 374? up to the required location (or I9839J is 

used for the "^trix and its ■nanipul^ti-n . The n-trix is stored in column sequence ->nd uses memory from 03749 upward. 

The following diagram indicates th? T.T-ner in which n^trix element": are stcredt 



00 



3760 






S+IO(2m+3) 


10/Cj 


S+I0((n+!) (m+2>-a> 




s 


Functional 


S+|0(2m+4) 




S+|0(tn+l)(m+2)» 






S+ !0(ro+3> 


S+IO(2m+5) 




S+|0((n+l)(m+2)+n 


il umn 








ncre^sing j 




3 

c 

1 

37^0+10(tiM-l) 


< 

Q 

S+ ICm 


S+IO(2m+2) 


o» 
c 

J 

S+IO<3m+4) 


j. S+|0(< J+1 )(m+2)*l ) 


S+|0((n+2)(m+2)-2) 



S Is computed for ench problem and Is equal to 03760 + l0(m+3) # All addresses are Maid addresses. 

2. PROGRAM ORDER (Card numbers are tn columns 78-80 of program decks) 

a. DATA LOADER 000-053 

b. COST CHANGER 100-107 
C. RHS CHJ1GER 200-260 



d. SOLUTION 300-916 (Not continuous) 

(1) Floating Point Subroutines 300-328 

(2) Shadow Price 400-423 
<M Duet 500-526 
(4) Simplex 600-628 
(-5K Final Basis 700-720 

(6) Fix and f rint 750-768 

(7) Non 3 a sis 600-821 

(8) Vatrlx Punch 900-916 



00-80 listings, v/lth remarks, of the machine language prpgr-m decks follow. These csrds are in a format similar to 
the condensed output. SpeciflcaMyt 

Columns Contents 

1-60 Instructions or constants 

61 * (record nark) if full 60 columns are us?d. Other wise the record mark is at end 
of the data being punched. 

62 I 

63-64 St Loader Code 

65-69 Address of where tf)s first character on the card is to be placed in memory. 

70-74 Address of wh;re the taSt character will f*»U In memory 

75 Flag 

76-77 Wot used 

78-80 Ca,rd sequence numbers 

TMs llst5ng was made on a 407 Accounting Machine. Since the 407 does not have flags, Ragged digits will appear as 

=tph*bet1caf characters in the J-R range. A flagged zero wi I I appear as a - (minus) or » zero. (The 407 was a 

bit recalcitrant.) A record mark prints out as a 2. 



4. in order to allow this program to run on a machine with greater than 20K memory, certain changes have to be made. Any 
digit In the attached list that is marked with an * must be changed. In each cose this digit will be a I (one) Oi 
a J (flagged one). This 1 (or flagged 1) must be changed to a 3 (or flagged 3) for 40 K memory, or a 5 (or 
flagged 5) for 60K. I.e., add two to the digit for •ach additional 20K of storage added. 



£>ITA LoAbc-fL. dvi Mos 00$ - 053 

3600072005003400201009004400012002752600099002742900011000002400090002*900000000 F.*»+ C*.*4 Data LoAOtO^ 
26000990026431000000020026001 1400274290000000011490001200009000000000000 - 1 
440000000000160231300001370231500500250232500400310004802314Z0010040200462000002 
190374200000360231400500250231700400^1037430231421005320374520010046200522000003 
1600064037901503746000002600079037461 103746000K0260006903746Z0010052200582000004 
230232100400310374302318210062803745160007400000250374502322Z0010036200642000005 
4300734037452 100684000643602314005003100001023141 100684000P020010064200702000006 
450O666023854800O000OO00491984O0OO00Z000O0O00O0000O - 10070200738000007 

3702315OO50025023510O4O0250239100400250243 100400250247 100400Z00 10073400794000008 
430095002316140231300004470087801200160231300001 160092502316Z0010079400854000009 
160094902314490073400000 11023130000 11 100925000401 10094900040200100854009140000 10 
4300938023 1649008060000031023 1402314250374002 3152503 741 023 17Z00 1009 14009740000 11 
250374202319320231800000 1402319000P04601 03401 300320374200000Z001009 740 10340000 12 
2503743023232503744023252503745023273202326000001402327000POZ0010 10340 10940000 13 
4601 11801300320374500000330231800000330232600000320374000000Z001010940 1154000014 
1403742000004602134012003203743000002303745000692 10009900069ZOO 10 11540 12 140000 15 
3300099000002 10009803742 2 1000980006349022 74000002602 12000099Z0010 12 140 12740000 16 
440147403745 160231 100000250041602331 11 013040000 11 101 30900002ZOO 10 12740 13340000 17 
110231 100001 140231 1000J0470129601 2003202348000001402 349000P0Z0010 133401 3940000 18 
46014180130032O042500O00330235000OO0160 1304004161601 3090233 1Z00 10 1394014540000 19 
32O0416OOOOO49O21140O000Z000000 -10145401478000020 
440 149403742490 1286000O0Z00O0O0 -10147401498000021 
2502351004003103726023303203726000001403727000P0460165401300ZOO 1014940155400002 2 



The data load routine rtid$ 
and stores parameters for 

the problem. It then pro- 
ceeds to load and store the 
matrix In the proper format. 
In the case of fixed point 
input, storage locations are 
computed from the row/column 
designations, and the Input 
elements are converted to 
floating point notations as 
they are read and stored. 



1403727000K0460163401200430165403727330166500000320372800000Z001015540 1614000023 

310372603728490 166600000ZOOOOOO -10161401638000024 

320166500000490160200000ZOOOOOO -10163401658000025 

330166500000 16004170O0N0 160172003727 1601 70803 7283201 70200000Z0010 16540 1714000026 

1401714*00003460 17820 1200110172000002 11 01708000021 10041700001Z0010 171401774000027 

49017O2OO0O0Z0000000O0OO00 -10177401786000028 

26018290170826018240 1720 12018240000 131018 1801 8184301 88603727Z00101782018420000,>9 

31037260372812004170000 1450183003727490080600000Z000000000000001018420 1890000030 

160146500009160194503727160196000418160196503 72745019540000020010188601946000031 

49O20220OOOO20OQOO0O0O0O00 -10194601958000032 

250000000000 110194500002 110196000001 110 196500002 120146500001Z00 10195402014000033 

490193400000ZOOOOOOOOOOOO - 10201402026000034 

260206401960 120146500001460209001200150205800000110206400001Z0010202202082000035 

4902034OOO0OZOO0OO0O0O0O00 -10208202094000036 

4402 1140 16653200425000002602 11400425490080600000Z0000000000000010209002138000037 

2602212000641202212000J9230006900073320009500000210009900069Z00102 13402 194000038 

2 100099000693202206000001 1022 12 000J01 200099000 J0470220601200ZOO 10219402254000039 

►480000000000491984000000ZOOOOOO ^102254022 78000040 

3303743000003200095 0O0OO49012620O00OZ0O00O0000000O0 -10227402310000041 

OOOOZOOOOOOOOOOOOOOO 1010231002314000042 

480000000000ZOOOOOOOOOOOO - 1023 1402326000043 

* * * *** * * 

36 1996500500441994400039 26 198670003825 1985 119876261990600033Z001J9840J9900000044 

3U990ot996526199300003825t9924f985l49?984000000Z000000000000001J9900J9948000045 



261996200038491995600000ZOOOOOO 



- 1J9944 J9968000046 
-10009600115000047 



L600000005004900000ZOOOOOOOOOOO 
36001000050036001720050036002440050036P031600500360000000500000000000000000 48 

000000000000102030400020406080003060902100408021610050015102006021814200Z0000049 

704112820080614223009081726300000000005060708090012141618151811242720242Z0000050 

822363520353045403632484455324946536048465462754453627180123456789123456Z0000051 

7890234567890J34567890jK4567890JRL567890^LM67890JKLMN7890JKLMNO890JKLHNZ0000052 

41000000000049004020P90JKLHNOPOZOOOOL100388000X9M90000000000M9000360000000000053 



iav£ D*TM LoAO«*- 



cost cH^)6-e/i- 



Qtvl Mov too- 



480000000O004602158O02003402511005004902182034000O0001023602Z101021 1402174000100 
51 100100320251 10000014025 14000004619840012004902254002003400Z1010217402234000101 



f.«.+ <r a v 4 co%r- cm vc~eji~ 



000001023803737001002602 193000792602296000641 102296000042400Z1010223402294000102 
000025144602448012001 102296000J01202193000J04702290012002602Z10 10229402354000103 
193000742602392022962102392000692400000025144602468012002102Z1010235402414000104 
392000691202193000J04702386012004826024980229649024800260249Z1010241402474000105 
8O239211O2498OOOO626OOOO0O25204902126ZZOO0O00OO0OOOOOO0O000 101024 7402512000106 
000000000000 -2114000107 ^i-r- CtU 



Cost afam-mL, 



The cost change routine rtads ID/cost elements, starches th# matrix for a corresponding 
identification number, and stores the new cost coefficient in the proper location. 



e-2 



1600536005384900434016005360055826011970000026011 1200445 1601Z1010040200462000200 
124000001 101 1240000 146004940 14004901 128004692601 17600000260 U 10 10046200522000201 
1870182749000000220 1 I860 1 1764900570Q260 1 18601 176 150084700001Z10 10052200582000202 
16008570118726010920118922011780118947007020110026011760119721010058200642000203 
260 119701186260 11860 11 76210109201 189260 11 7801 1B9320117800000Z10 10064200702000204 
1501 198000001401 1780000Q470 109401 300440077<jOl 197320 119800000Z 10 1007020076200020 5 
330 1197000002 10085701 178440082201 186330 118600000150084700002Z 10 10076200822000206 
150117800000150118900000210119800000250109301 1983301 19800000Z 10 10082200882000207 
430098601 189430103401 190460096601 200110 119001 1912501 1980055 7Z 10 1008820O9420OO20 8 
,120109200001460089401300260119 7018264901 10602601 1980 11 97330 1Z 10 1009420 1002000209 
190000001101092000014601 1260 1400150 118 9000001 10 1 198000054300Z 10 10 1 00201 0620002 1 
9860 1189440109401093320 119700000260 11 890 109226000000 11974900Z 10 10 10620 11 220002 11 
000048260115100469260051700000490050600000000000000000 Z 1 1 1 1220 1 1 820002 1 2 
000O00O0OO00000O02Z000000000000 1010118201200000213 
160053601256 150 175300002 1601 112000993201 1900000049004580430 1Z 10 10 1200012600002 14 
^ 27601 19069009660430 12900 1169482201 1890 1168250 1325011 763301 17Z 10 101 2600 13200002 15 
I 6000001501 16800000 160134700000 11013470000 12 10 11 76011 764301 39Z10 10 1320013800002 16 
40 1 1684901 3500160 14410 1837250 14280 11 69 1101441 00000230 1175000Z 10 10 13800 14400002 17 
0026014770144 11 201477000042601 166000002201 166000933201 162000Z 10 10 14400 15000002 18 
002301 1760 11652601 1860 19372201 18600096320 11 780000023011 8601 1Z 10 1015000 15600002 19 
863300082000001 50008 10000J2 10009001 186 2601 18600090230 11 8601 1Z 10 10 15600 1620000220 
652 10009500095 120 13470000 14701622012002601 18600094440169401 3Z1010 162001680000221 
253201 186000002301 18601 1971 10 11890OONO250 1093000992601 09201 1Z 10 10 16800 1740000222 



The requirement change routine reeds 
variable Identification numbers end 
changes for the associated requirement 
element. The chaeeje it converted to 
floating point notations, the metrlx 
Is searched for the Identification 
number, and the requirement vector Is 
updated. 



892601 198000924601774013004900966047008820 14004901 1260000000Z 10 10174001800000223 
0100000000000 - J90692R080J74077P568J601280405J4Z 1010 180001860000224 

8249H490J38013M759J29099M164J21267L674J14332L266J08147K922J0Z1010 186001920000225 
2597K630JOO0O0OOO01 60053601970 150175300003490 12240430199001 1Z 1010192001980000226 
90490096.60430201001 169490096603201 16900000 120 1168000N1230119Z1010198002040000227 
701176210118901 16849O1718ZZO0O000OO 1010204002066000228 
480000000000460221800300370348500500250229203485250229303487Z001021 1402174000229 
250229403489250229503491250229600400490224200000Z0000000000000010217402222000230 
340O00O0O10236O2292001003202292000O0 14022950000046034640 1200Z0010221802278000231 
470232200300340000000 108 370372700 10049024 1400000ZOOOOOOOOOOOOOO 10227802326000232 
320350000000 140350 1000P0460239001 30043023900350 1250352 100400ZOO 10232202382000233 
49024O20OOOOZOO0O00000O00 - 10238202394000234 

250351900400310372603500320372600000 1403727000P046025500130020010239002450000235 
1403727000K046025300120043025500372733021370000032037280000020010245002S10000236 
3103726037284902562 0OOO0Z0O00OO -10251002534000237 
3202 13700000490249800000Z 000000 -10253002554000238 
3302 13700000 16021 17000N0 160261603727 160260403728320259800000200 10255002610000239 
1402610000034602676012001 10261600002 1102604000021102 11 70000 1Z001026 1002670000240 
490259800000Z000000000000 -10267002682000241 
260272502604260272002616120272000001310271402714430278203727Z0010267802738000242 
310372603728120211700001450272603727490212600000Z0000000000000010273802786000243 
16022530000916028410372716028560211816028610372745028500000020010278202842000244 
49029180OO00Z0000000OOO00 - 10284202854000245 



f-1 



2502850028501102841000021 10285600001 110286100002120225300001200102850029100002*6 
4902830000002000000000000 -10291002922000247 
26029600285612022530000146029860120015029540000011029600000120010291802978000248 
4902930000002000000000000 -10297802990000249 
44030100213732021250000026022530007926030520006411030520000420010298603046000250 
240304602295460319201200 1103052O00J01202253000J047030460 1200200 10304603106000251 
260225300074260314803052210314600069240314202295460327201200200 10310603166000252 
1202253000J0470313001200480000000000200000000000000 -10316603202000253 
26032510305221032510006911032510000616004690326316004450324020010319263252000254 
4900422021254902126000002000000 -10325203276000255 
2603355031481103355000J6260339100064210339100069260225300079Z0010327203332000256 
160046903367^6011970334449019380212516004690340316004450338020010333203392000257 
4900422000991103355000J01103391000J01202253000J0460333201300Z0010339203452000258 
49021260O000480OOOO000OO49198400000OZ000000O0000000 -10345203488000259 

000000000000 -2114000260 < M f W&-€JL~ 



160053600538490043401600536005582601197000002601112004451601Z1010040200462000300 aitt fSu^^ S«kw*r*4, 

12400000 1101 1240000146004940 14004901 12800469260 1176000002601Z10 10046200522000301 
18701827490000002201 1860 1176490057002601 18601 17615008470000121010052200582000302 
1600857011872601092011892201178011894700702011002601 17601 197Z1010058200642000303 
*60119701186260118601176210109201 18926011 7801 1893201 17800000Z 1010064200702000304 
150119800000140 117800000470109401300440077401 1973201 19800000Z1010070200762000305 
3301197000002100857011784400822011863301 18600000150084700002Z 1010076200822000306 
1501 17800000 1501189000002101 19800000250 109301 1983301 19800000Z10 10082200882000307 
430098601189430103401190460096601200310119001 1912501 19800557Z1010088200942000308 
1201092000014600894013002601 19701 826490 11060260 119801 197330 1Z 1010094201002000309 
190000001 101092000014601 1260 14001501 189000001 10 1198000054300Z10 10100201062000310 
98601 1894401094010933201 197000002601 18901092260000001 1974900Z10 10106201 122000311 
0OOO4826OH51OO4692600517OO0O04900506000OOO000000OO00O00O ZlO 10112201182000312 
000O000OOO00O0O00ZZ00O00OO0O0 1010118201200000313 
160053601256 1501753000021601 112000993201 1900000049004580430 12 10 1012000 1260000314 
27601 1904900966043012900 1169482201 1890 1168250 132501 176330117Z1010126001320000315 
6000001501 16800000 1601347000001 10 13470000 12 10 11 7601 176430139Z 1010 132001380000316 
40116849013500160144101837250142801 1691 101441 000002301 175000Z1010 138001440000317 
002601477014411201477000042601166000002201166000933201 1620OOZ10 10 144001500000318 
002301 17601 1652601 186019372201 186000963201 178000002301 18601 U 10 10 150001560000319 
863300082000001 5000810000J210009001 1862601 186000902301 18601 1Z10 10156001620000320 
6f>2H000950009512013470000 14701622012002601 186000944401694013Z 1010162001680000321 
253201186000002301 18601 1971 101 189000N0250l093000992601092011ZlQiClv>8e? 1740000322 



The shadow price calculator converts 
all cost coefficients to floating 
point notation and evaluates Zj -cj 
for each non-basis variable. 



892601 198000924601774013004900966047008820140049011260000000Z10 10174001800000323 
000000000000- - J90692R080J74O77P568J60128O405J4Z1010180001860000324 
8249H490J38013M759J29099M164J21267L674J14332L266J08147K922J0Z1 010 18600 1920000325 
2597K630J00000000016005360 19701 501 75300003490 12240430199001 1Z 10 10 19200 1980000326 
9049009660430201001 169490096603201 16900000120 1168000N1 2301 19Z10 10 198002040000327 

70117621011890H684901718ZZOOOOOOOO 1010204002066000328 Ui+ Cnl - Floa-k™, 

1602472037802602173000642602253000791102173000J0260375800000Z1010211402174000400 rW Ca^ - Shwiou, P*nc 
320375300000220375801822470223001 200260376001 82649024660 1603Z 10 1021 7402234000401 
752000N33303753000001602337037531602429018181602436037592503Z10 10223402294000402 
75003758330375800000250375903 1554302394000001 1023370000 11203Z 10 10229402354000403 
75200001110242900001 1202436000014902326026024170233731037530Z10 10235402414000404 
000026037600 18183303759000004402466037 50320376000000 26000000Z 10 1024 14024 74000405 
37601 102472000J01202253000J04702 1500 12002602472000642 1024720Z 10 10247402534000406 
00692102472000691 1021 73000J026021610006926024890006926022530Z 10 1025340259400040 7 
00741602448026481602455000M1 2602520026474902 150M902668320376Z 10 10259402654000408 
000000490246602602659000 7426027 10000642 1027 10000692600000018Z 10 10265402714000409 
26260287 1027 10260284702710 1102847000J01602799037802602253000Z 10 1027 14027740004 10 
792603760018262 1037600000046028840 1200 1600469028472601 197037Z 10 102774028340004 11 
60490 193800000 16004690288 3 160044500000490042200099 11 02 799000Z 10 102834028940004 12 

J01102847000J01202253000J047027760 1200 1102847000K02 10287 1000Z 10 102894029540004 13 

* * 

691 202659000 J046027520130036 1996400500340000000 10239 1996500 1Z 10 1029540301 40004 14 

00390004900 100340000000 102340000000 102 16026590000436?9964005Z 10 103014030740004 15 
^ 0047031 12001 0039f996500 1003400000001 08 1202659000014703064012Z 10 103074031 34000416 



OO340OO0O0010249T9840ZZ0O000OO0 
0000000000000 
4341624500OZ00OO000000OO0 
49634559OO55560ZOOOOO0O000O0 
4464554 36349565541530ZO0OOO0000 
0000654159005664630Z 10000000000 



1010313403156000417 
-2114000418 
419 
420 
421 
422 



OOO06541590049550Z0OOOO000O0OO0 423 Cl^ - ^uiJo*, A» 

2600023018262603440000792602185000642102185000691102185000JOZ1010211402174000500 - Dv«.l 



240002300000470222201 100260222 1021852600023000001203440000JOZ1010217402234000501 
-^47021620120044?9840000232600023022051 6035600000026034 1600074Z 10 10223402294000502 
260237202221260241300064210241300069210237200069260242502372Z10102294023S4000503 
2102413000692400000018264602486013001600469024252601 19700000Z 1010235402414000504 
490120OOO00024O009900023470248601100260002300099260352402372Z1010241402474000505 
260356002413 1203416000J0470233001200140356000000470254601200Z 10 10247402534000506 
4M110000000260268403524260270803560260344000079260266503524Z1010253402594000507 
160269603770260278003524260270103560 1600469026652601 19702545Z 1010259402654000508 
490120000000260376000099260000001826260377000000260000001826Z1010265402714000509 
1102696000J01102701000J01102708000J01203440000J0460269001300Z10 102714027740005 10 
26000000254526028570222 1260286402221260341600074160046902857Z10102774028340005 11 
2601 197037604901938000002600000000992102857000692 102864000692101028 34028940005 12 
1203416000J0460282201 300 160298403770260302500064210302500069Z 10 1028940295400051 3 
160309703770260344000079220377001826460321801 20026030850222 1Z 10102954030140005 14 
160312100000260341600074260314003121 120314000008 160046903097Z10 10301403074000515 



The dual algorithm computes a feasible 
solution for the problem. Its operation 
may be monitored allowing the typing of 
the value of the functional and the 
Identifications of the variables being' 
entered Into and removed from the basis. 
It should be noted In this regard that 
the objectives of feasibility and 
optimal Ity may not be compatible at 
this point, and therefor*, the functional 
way actually decrease from one iteration 
to the next while the dual algorithm is 
In control 



During the reduction of the matrix of 
coefficients a test is made to determine 
whether the resulting quantity has a 
magnitude less than an amout cetled 
'essential zero.' Any element falling 
into this category is set to z ero In 
order to avoid computation with numbers 
which are not significant. 'Essential 
zero' has been set to 10" 8 . Its value 
(the corresponding floating point 
exponent) is located in memory addresses 



2601 19700000490193803770 160046903 13 3 16004450000 0490040200099Z 10 10307403 1340005 16 
1400000000M34703598013002 103085000692103140000692103 121000692 10 103 13403 19400051 7 

03144 and 03145 (Card Mo. 517) and 

1203416000J04603062013001102984000J01103097000J01103025000JOZ1010319403254000518 may be changed If desired* If this 

change Is made, It Is Important to 

1203440OOOJ04602978013O026O3373O222126033610356.0120336100OJ0Z10103254033140O0519 note that the high order digit must 

carry a flag. 

260336803361220337300069260338003373260002300000260000000000Z 1010331403374000520 
260000000023 1103631 0000147021 1400100340000000102380362900100Z 10 10337403434000521 
340000000108260357703373260354103361260350500064210350500069Z 10 10343403494000522 
260374600000380373700100340000000 108260374600000380373700 100Z 1010349403554000523 
34000000010826037460000038037370010049021 14026036160312 12600Z 10 10355403614000524 
OOOO182649O3158OO0ZZ0O0OOO000O0 1010361403633000525 

000000000000- -2114000526 JUs4 Ox^ri - Out*, i 

250366803747260002301826260354800074260219700064210219700069Z1010211402174000600 C^U - S.^pUy 

2102 19700069240002300000470223401100260223 302 197260002 300000Z 1010217402234000601 
->2l02197000691203548000j047021860l20044r984000023260002303522Z 10 10223402294000602 

The simplex algorithm starts with any 
160348100000260363200079260238002233260243302233260242100064Z1010229402354000603 feasible solution anc! computes the 

optimal feasible solution. Its operation 
21024210006949025420240000001826470251801 1001600469024332601Z 1010235402414000604 may be monitored by allowing Iteration 

typeout as in the case of the dual 
197000004901200000002400099000234602518011004602586012002602Z 1010241402474000605 algorithm. 

493023802603746000002600023000992603481024211203632000J04602Z 1010247402534000606 The discussion of 'essential zero 1 as 

applied to frte dual algorithm is equally 
630012001102380000J01102421000J01102433000J04902374026026090Z1010253402594000607 valid for the simplex algorithm.; In this 

case, the test value is located in 
238024037460000047024700130049025180140348100000470266601200Z1010259402654000608 memory addresses 03264 and 03265 (Card 

48N110O0O0O02602804024932602828022332603632OO07926O278502493Z 1010265402714000609 

16028 1603770260290002493260282 102233 160046902785260 119702665Z 10 1027 14027740006 10 

490120000000260 37600009926000000182626037700000026000000 1826Z 10102774028340006 11 



: Jo. 619) 



In the case of a tie In the quotient, 
normal ly used as the criterion .for 
choosing the variable to leave th*. 
basis on any iteration, the choice 
Is resolved by using the largest 
denominator as a secondary criterion. 



1102816000J01 102821000 JO 1102828000 JO 1203632000J0460281001300Z 10 10283402894000612 
260000002665260297703481260298403481260354800074160046902977Z1010289402954000613 
2601 197037604901938000002600000000992 102977000692 10298400069Z 10 102954030 14000614 
1203548000J0460294201300 1603 104037702603145000642 103 14500069Z 10 10301403074000615 
160321703770260363200079220377001826460333801200260320503481Z1010307403134000616 
160324100000260354800074260326003241 120326000008 16004690321 72 10 1031340319400061 7 
260U9700000490193803770160046903253160044500000490040200099Z10 10319403254000618 
1400000000M3470370601300210320500069210326000069210324100069210 10325403314000619 
120354>SOOOj04603i82013001103l04000JOll03217000J01103145000JOZ10l033l403374000620 
1203632000J04603098013002603469022331203469000J0260347603469Z10 10337403434000621 
22034810006926034880348 1260002300000260000000000260000000023Z10 10343403494000622 
11036470000l4?b212600l60340000000102380366500 100340000000 lOffZ'IO 10349403554000623 
260368503481260364903469260361300064210361300069260374600000210 10355403614000624 
3803737OO1O034O0OO00010826O37460O00O38037370010034O0O0000108210 10361403674000625 
. 260374600000380373700100490212602603724032412600000018264903Z 10 10367403734000626 
' [7* 2780OOO0OO000ZZ000O00000O0 1010373403748000627 

OOOOOOOOOOOOO -2114000628 U±4- C u *4 S'VKpU? 

->491964000300340010000102260264500064210264500069361996400500Z1010256202622000700 CaA V'»u\ 

391996500100261998400000340001000108L8f9975001p0340000000l02Z 10 10262202682000701 
""^ 34000000010216035250000636*99640050039*9965001001203525000012 10 10266202742000702 
470270601200340000000102260331 3000641 103313000J0260336102645Z1010274202802000703 
1103361000J02603005033612103005000792602993030051102993000JOZ1010280202862000704 
2602957030051 102957000K0260353700079260075002658260074002585Z1010286202922000705 



16007250000016007150000016030360000026030890295726031010295721010292202982000706 
16030770OO00160314900000260320503005260352500074220O00003698Z 10 10298203042000707 
4603206012001600469030892601 19700000490 1200000004403 15800000Z10 10304203102000708 
240009900750470320601 10O26O075000099260073O00O004903206024002 10 10310203 162000709 
099007404603206013002600740000992600720000002 1030360006921032 10 10316203222000710 
077000692103089000692103 101000692103149000692 1032050006912* 3Zi0l032^26328200071 1 
52500OJ04703030012002619982000003400000O01023819973001003201Z10 10328203342000712 
955000002701948000002600800007143400000001083800797001002603Z10 1033420340200071 3 
421033132600760000003200755000004403478007603300760000003203Z10l0340'2b34620o67i4 
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RACD 
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00048 
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INPUT-3»AREAZ 







00510 


21 


00532 
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TD 
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BD 
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21 
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ARE AZ 
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-0001 


02314 
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TR 
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11 
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AW 
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45 
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02385 
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48 
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00734 
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43 
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BD 
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ALOAD 


CM 
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02313 
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CF 
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32 02346 
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01394 
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01418 
01430 
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00069 
00000 
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00099 
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02331 
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000-2 
000-1 
OOOJ0 
01200 
00000 
OOOPO 
01300 
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BE 
SF 
M 
A 

CF 
A 
A 
B 

TF 
BNF 

INITI TFM 
TRID TD 
AM 
AM 
AM 
CM 
BNE 
SF 
CM 
BNL 
SF 

BCLFL CF 
TFM 



OUTZ 
IMPUT-32 

IMPUT-1»MAM022 

99.MAND2Z 

99Z 

98* INPUT-4Z 



98fSA0DR-lZ 
XCLRZ 

ST0RE+6»99Z 
ROWMI 1»INPUT-1Z 
I .0.10Z 

EXP-1»AREA+17Z 

TRID+6t 1.10Z 

TRIO+11.2.10Z 

1*1* 10Z 

I .10.10Z 

TRIOZ 

AREA+34Z 



AREA+35t70.10Z 

BCLFLZ 

EXP+8Z 

AREA+36Z 

TRID+6tEXP-lZ 
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01454 32 
01466 49 
01474 
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01494 
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01506 31 
01518 32 
01530 14 
01542 46 
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00416 00000 
02114 00000 

01494 03742 
01286 00000 



02351 
03726 
03726 
03727 
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02330 
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OOOPO 
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03727 
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01578 43 
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01634 
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03728 
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TFM 
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B 

DORG 
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TR 
SF 
CM 
BNL 
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BE 
BD 
CF 
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TR 
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SF 

B 
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33 01665 
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32 01702 
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11 00417 
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TFM 
TFM 
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DEC CM 
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AM 
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B 
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SHFTDC TF 
TF 
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TR 
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TR 
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FLDDEFZ 

*-3Z 
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16 


01945 


-3 U i 




TFM 
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45 
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25 
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11 
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AM 
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11 
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26 
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TF 
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CTRtlt 10Z 
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AM 
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44 
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01665 
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32 
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00000 




SF 
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26 


02114 


00425 


STORE 


TF 
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49 
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26 
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TF 


FLG+6«SADDRZ 
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02212 


000 J9 




SH 
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23 


00069 


00073 




M 
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02170 


32 


00095 


00000 




SF 
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, 02182 
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00069 




A 
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02194 


2 1 


00099 


00069 




A 


99 »MAND2Z 


02206 


32 
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ouooo 


FLG 


SF 


*Z 
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} 1 


02212 
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AM 


FLG+6tl0.lOZ 






00099 
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5M 
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02242 


47 
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01200 
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FLGZ 


02254 


48 
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LSTARTZ 












DORG 


#-32 


02274 


33 




03743 


00000 


XCLR 


CF 


1NPUT-3Z 


02286 


32 


00095 


00000 




SF 


95Z 


02298 


49 


01262 


00000 




B 


INITI-24Z 


03760 
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DS 


t3760Z 


00079 
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DS 


• 79Z 


00074 








N 


OS 


t74Z 


00069 
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DS 


• 69Z 
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SAOOR 


DS 


• 64Z 
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02311 
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19944 
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25 19924 
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•0033 
19965 
•0038 
19851 
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26 19962 #0038 
49 19956 00000 



INPUT 


DS 


♦3746Z 


CASE 


DS 


♦ 58Z 


SIGN 


DS 


tCLRFLG+HZ 


CTR 


DS 


•FLSET+11Z 
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zooo 


0000 
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TF 

TD 

TF 

TR 

TF 

TD 

B 

DORG 

TF 
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LSTART+11»*Z 

LTR+6.LAREA+68Z 
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LSTARTZ 
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02114 
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2114Z 
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H 


Z 





02126 46 


02218 


00300 


GO 
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READTYZ 





02130 37 


03403 


00500 


READCD 


RACD 


RAREA+1Z 
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02190 25 


02292 


03485 




TP 
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02162 25 


02293 


03487 




TD 


IDIM-2.RAREA+3Z 





02174 25 


02294 


03489 




TO 
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